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FINAL REPORT 

HIGH TEMPERATURE FURNACE MODELING 
AND PERFORMANCE VERIFICATIONS 


1.0 INTRODUCTION 

Under NASA Grant NAG8-708 we performed analytical, 
numerical, and experimental studies on two classes of high 
temperature materials processing sources for their potential use 
as directional solidification furnaces. The research 
concentrated on a commercially available high temperature furnace 
using a zirconia ceramic tube as the heating element and an Arc 
Furnace under design at NASA's MSFC, based on a ST International 
tube welder. The zirconia furnace was delivered, assembled, and 
analyzed for its potential to function as a directional 
solidification furnace. The work, on the Arc Furnace was 
initially stalled due to the unavailability of the NASA 
prototype. A proposal was written and funded to purchase an 
additional arc welder to alleviate this problem. The ST 
International weld head and power supply were received in 
September 1988 and used to experimentally verify a portion of the 
numerical results generated by a thermal model using SINDA. 


2.0 OBJECTIVES 

There were five major objectives researched during the 
grant. The first objective was to assemble the zirconia furnace 
and construct components needed to successfully perform a series 
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of experiments. The second objective was to evaluate the 
zirconia furnace performance as-delivered from the manufacturer 
to evaluate its potential as a directional solidification furnace 
element. The third objective to establish a data base on 
materials used in furnace construction, with particular emphasis 
on emissivities, transmissivities and absorpt ivities as functions 
of wavelength and temperature. These data were used to develop a 
one-dimensional and two-dimensional spectral radiation heat 
transfer model for comparison with standard effective radiation 
heat transfer modeling techniques. These models were used to 
predict wall and crucible temperatures and may serve as an 
effective method to apply open loop control techniques to 
advanced high temperature resistance and arc furnaces. The 
fourth objective addressed the development of a SINDA model for 
the MSFC Arc Furnace and was used to begin preliminary numerical 
designs of sample holders and to estimate cooling media 
temperatures for the steady-state operation of the furnace. 
Finally, the fifth objective addressed the initial performance 
evaluation of the Arc furnace and associated equipment for 
directional solidification. 

These objectives were met by the research. In addition, 
preliminary research was conducted to determine if a zirconia 
heating element could be developed that would function longer 
than the two commercially prepared elements purchased with the 
ARTCOR Model 460-15 Bench Top Tubular Laboratory Research 
Furnace. A zirconia heating element was prepared and its 
performance reported. 
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Objective 1 . 

One of the focuses of this research project was the 
evaluation and performance verification of existing furnace 
technologies. Two furnace types are being considered. They are 
the ARTCOR Model 460-15 Bench Top Tubular Laboratory Research 
Furnace and an Arc Furnace being constructed at NASA's MSFC which 
uses an ST International Orbital Tube Welding System as a heat 
source. This final report's first three sections concentrate on 
the ARTCOR Research Furnace. 

Objective 2. 

The second section contains an initial performance analysis 
of the zirconia furnace. The performance of the zirconia furnace 
was tested at different levels of applied power to the main 
zirconia heating element, namely 15%, 20%, 25% and 30%. The axial 
temperature profiles in the furnace were measured along the 
centerline for applied powers of 20% and 30%. Step change 
increase from 20% to 30% applied power and a linear ramp rate of 
1% per minute from 0% to 15% applied power were studied. 
Limitation of 50% applied power were corrected with repair of the 
units controller. Power to higher levels was then applied. 

Objective 3 . 

Experimental data on various ceramics and other materials 
needed for eventual design, construction and modeling of advanced 
high temperature electrical furnaces were collected from various 
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]_ it.6rsit.ur6 references. Spectral radiation data were collected on 
alumina, zirconia, and titania. The data were assembled to 
evaluate the material's emissivity, transmissivity, and 
absorptivity trends at various wavelengths and operating 
temperatures. Data for zirconia were used to estimate the 
effective emissivity for use with the high temperature optical 
pyrometer. Finally, the thermal and frequency characteristics of 
these data were used to formulate a one - dimensional and 
two-dimensional spectral model for comparison with effective 
radiation modeling approaches. FORTRAN computer programs were 
developed to solve these highly nonlinear problem. These efforts 
were coordinated with those of the UAH Center for Microgravity 
Studies to prevent duplication and to accelerate development of 
advance high temperature material processing technologies. 

Objective 4 . 

Numerical Analysis of the MSFC Arc Furnace. An Arc Furnace 
for use on the KC-135 experimental aircraft was first proposed by 
Poorman (1986). The furnace uses an S.T. International orbital 
tube welding head and power supply to produce localized and 
extremely high temperatures. Several preliminary tests were 
conducted by Poorman on samples of Al, Cu and W using the 
configuration shown in Figure 1. We numerically characterized the 
processes and design parameters involved in the development of 
the Arc Furnace for use as a directional solidification furnace 
at temperatures exceeding those currently available, i.e. 1600 C. 
This portion of the research was directed towards developing 
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scenarios for sample support, external furnace shielding, 
translation mechanisms, temperature limitations, and extensions 
of numerical characterizations to nonmetallic materials. 

A numerical model developed for the Arc Furnace and 
described in following sections and was applied to a sample 
holder with physical properties of stainless steel, aluminum and 
copper. The goal of this numerical study was to estimate the 
performance of these holders under various isothermal external 
boundaries. The external boundary was held constant at 
temperatures consistent with air, water, and oil cooling (10, 38, 

50, 100, 150 and 200°C) . Axial and radial profiles were 
determined for 75 mm long samples of tungsten held at each end by 
25 mm long aluminum, copper, and stainless steel holders. 

Objective 5 . 

Performance Verification of the Arc Furnace. The S.T. 
International weld head and power supply were received in 
September 1988, and after several repairs, was placed in service 
in February 1989. Tests were conducted on 3” stainless steel and 
aluminum tubing while using a specially designed water cooling 
system. The water cooling system keeps the weld head cool during 
high temperature operation associated with welding large diameter 
specimens at high currents. Details of this system are discussed 
in Section 3.5. 
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3.0 RESULTS AND DISCUSSION 

This research supported a graduate student that prepared a 
thesis on the numerical modeling and experimental testing of the 
ARTCOR Model 460-15 Bench Top Tubular Laboratory Research 
Furnace. This thesis documents the work on this portion of the 
research and is attached as Appendix 1. The following sections 
directs the reader to the sections covering the specific research 
on each topic. 

3.1 Objective 1. ARTCOR Furnace Assembly. 

The ARTCOR zirconia furnace was received from NASA and 
required substantial assembly. A stand was constructed to 
support the furnace, optical pyrometer, and ventilator. This 
stand was designed to protect the furnace from vibrations and was 
fitted with clamps to hold accessories needed during the 
experiments. Detailed procedures are given in Chapter 3 of 
Appendix 1 under the topic commercial furnace. In addition, a 
research version was developed which was used to research the 
furnace and element performance. Information on the research 
version of the furnace is also provided in Chapter 3 of Appendix 

1 . 

3 . 2 Objective 2 . Zirconia Furnace Performance Verifications . 

One requirement of this research was that the commercial 
zirconia furnace initially be evaluated in an as-received 
condition. No efforts were made to improve, modify, or alter in 
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any way its geometry, heat transfer characteristics, or control 
routines. It should be noted that free convection is used to 
cool the furnace casing. Also, since the furnace ends are open to 
the atmosphere, one recognizes that free convection and radiation 
losses will play major roles in the furnace's performance. 

Attempts to utilize this furnace for materials processing m a 
micro-gravity environment will require active cooling to correct 

this limitation. 

A major problem with the commercial zirconia furnace was in 
the area of control. The furnace had no feed-back mechanism for 
thermal control or even a direct readout for furnace temperature. 
These problems were outlined in our semiannual and annual 
progress reports on this research. Basically the controller was 
found to be ineffective for controlling the furnace element, such 
that scientific temperature data could not be collected. No 
technique was available with this equipment to set or control the 
furnace to a desired operating temperature. In actual operations 
with this furnace, direct measurement of the temperature within 
the furnace and manual adjustment of the percentage power applied 
by the controller to the furnace are required get the furnace 
operating properly. Accurate temperature control is required to 
develop this type furnace for eventual directional solidification 

experiments . 

The furnace support structure for the element did not 
produce dimensional stability that would facilitate consistent 
clearances around objects placed within the element. The 
dimensional limitation and measurement of the temperature 
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profiles were addressed in this research as part of the research 
furnace and reported in previous progress reports and reviewed in 
Chapter 3 of Appendix 1. In Appendix 1, some of the possible 
ways to correct the variations in element placement within the 
furnace are discussed. Some of these corrective techniques were 
incorporated in the research version discussed in the same 
chapter . 

Element failure due to crack development were experienced 
with the commercial supplied elements. These problems were 
reported in previous progress reports and are reviewed in 
Appendix 1, Chapter 3. A research element was developed and 
cured in such a way that extended the life time of the element. 
Procedures used to prepare this element and measure the power 
consumed by the furnace and thermal distributions are discussed 
in Chapter 3, Appendix 1. 

Detailed results on the above tasks can be found in our 
previous semiannual progress reports, the first annual report, 
and in Chapter 4, sections 4.2 of Appendix 1. The major results 
of this research show that the current configuration of the 
ARTCOR Model 460-15 Bench Top Tubular Laboratory Research Furnace 
lacks sufficient controls to provide scientific temperature data 
and dimensional stability to retail the element in a fixed 
location. Lack of element support on one end is an engineering 
problem that can be addressed with simple modifications as shown 
by using the research version prepared during this research. 
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The control problem is a larger problem that requires 
further research. The upgrade of the controller or the inclusion 
of a computer/controller combination that would permit feed back 
control based on temperature measurements would greatly improve 
research potential of this commercial furnace. 

3.3 Objective 3. Zirconia Furnace Model Development. 

In order to improve the performance of the furnace, an 
accurate model for the furnace has to be developed. Detailed 
models were developed for comparison with routine radiation heat 
transfer models which uses bulk radiation properties. A 
literature search was conducted on the spectral properties of 
furnace materials of construction namely alumina, zirconia, and 
titania. In addition, thermal properties of iron were also found 
and used in the numerical analysis. Radiation properties, which 
are functions of wavelength, were used to predict the radiation 
heat transfer of the crucible wall. At issue was the direct 
coupling of the sample to the furnace wall at specific 
frequencies for which the crucible wall is transparent. The 1-D 
and 2-D models take into account the absorption and conduction of 
the crucible wall to the iron sample. 

Development of all aspects of this portion of the work along 
with a review of other measurement techniques are provided in 
Chapter 1, Chapter 2 (formulation of the one- and two-dimensional 
spectral radiation model) and Chapter 4 section 4.1, where the 
numerical results may be found. 
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The major finding of the numerical studies is that spectral 
radiation transfer plays a major role in zirconia type furnaces, 
and that with proper selection of the radial dimension of 
components used within the furnace, thermal gradients and their 
resulting losses can be minimized. Two papers, one on the 
one -dimensional model and its results and another on the 
two-dimensional model are in preparation. The one-dimensional 
model was submitted for review and returned for correction. It 
is currently in rewrite. 

3 . 4 Objective 4 . Arc Furnace Model Development . 

A numerical model of the Arc Furnace was been developed to 
evaluate a conceptual design based on the work of Poorman. The 
studies were made using an idealized model with configuration 
shown in Figure 1, but with a sample rod only 20 mm long. Steady 
state idealized models of the Arc Furnace were numerically 
simulated using SINDA. This numerical model assumed symmetry; 
therefore, only half the configuration was modeled. 

Additionally, temperatures in the circumferential direction were 
assumed constant. The numerical model consisted of a 6.35 mm 
O.D. x 70 mm long sample of pure tungsten rod with 25 mm of rod 
inserted into each holder of 13 mm O.D. x 25 mm long stainless 
steel. This configuration provides only a 20 mm long tungsten 
sample rod. Since, symmetry was assumed, only 120 mm of the 
sample rod and one heat sink required numerical modeling. We 
limited our numerical analysis to a sample length of 70 mm, since 
this is our current experimental limit for preparing such a 
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The numerical and 


sample by powdered metal sintering techniques, 
experimental capabilities can easily be extended to a longer 

sample length. 

The model developed included an isothermal boundary imposed 
at the end of the sample rod. The isothermal temperature was 
maintained at the melting point of pure tungsten (3387°C) . The 
isothermal surface of the stainless steel holder was set at the 
temperatures consistent with air, water and oil cooling 
(10-200°C) . Additionally, it was assumed that there was zero 
contact resistance between the tungsten sample rod and the 
stainless steel heat sink. An idealized radiation boundary was 
imposed on the exposed surface of the sample. This surface was 
radiation coupled to a constant 45°C environment and was assumed 
to be a sufficient distance away from this environment to make 
the view factor equal 1. The emissivity of tungsten at 3315. 5°C, 
was assumed constant at 0.39 (Holman, 1981). To simplify the 
model, radiation exchange between the sample and sample holder 

was initially neglected. 

The model was developed on a 1 mm grid spacing which 
required 222 nodes, with additional nodes for radiation to the 
environment. The steady state solutions took approximately 3 
minutes of computer core time. The actual performance of the 
experimental furnace, processing a tungsten sample, should fall 
below this idealized case, due to the steady state assumption. 

The power consumption as a function of the temperature of 
the quench block is shown in Figure 2. Power consumption ranges 
from 1.76 to 1.70 KW at a quench block temperature of 10 and 
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Figure 2 . Steady-state power required for different quench block 
temperatures . 


- 13 - 



200°C, respectively. An increase in power of 0.32 watts is 
required for each degree centigrade when the quench block 
temperature is operated below 200°C. Values discussed here and 
in the figure are within the experimental range determined by 
Poorman (1986) for his 3.2 mm diameter tungsten sample. 

Centerline axial temperature profiles for tungsten samples 
calculated at various quench block temperatures are shown in 
Figure 3. An axial centerline gradient of 2650°C was calculated 
for the first 3mm of the sample. Figure 3 shows that this 
gradient is virtually constant for the first 2-3mm on either side 
of the arc location for the quench block temperatures imposed. 

The holder's isothermal boundary begins to influence the axial 
gradient in the sample at a distance of about 5 mm from the arc's 

location at 0 mm. 

Several isotherms for the model were included in the contour 
plots shown in Figure 4. These temperature contours range from 
the quench block temperature to 3387°C and are plotted at their 
respective positions in the modeled portion of the sample. The 
calculations indicate that a relatively small, and extremely 
localized melt zone exists. In the first millimeter along the 
axial direction from the melt interface the temperature drops 
187oC from the melting point of tungsten at 3387 C to 3200 C. 
Figures 4 a, b, c, and d show contour plots obtained for the Arc 
Furnace with isothermal boundaries of 50, 100, 150, and 200°C, 
respectively. In all four figures, the 500°C isotherm penetrates 
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Figure k. Contour plots obtained for the Arc Furnace with a 

tungsten sample and stainless steel holder. Boundary 
conditions on the external surface of the holder were 
controlled at a) 50°C, b) 100°C, c) 150°C, and d) 
200°C. 
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well into the holder region. This temperature is approaching the 
melting point for aluminum holders. Even at temperatures as low 
as 10 and 38°C see Figures 5, this penetration is still apparent. 

In all six figures one consistent trend can be seen. The 
isotherms in the samples from the center of the sample (0 mm in 
Figures 4 and 5) to about 10 mm are relatively constant no matter 
what cooling fluid temperature is used. This trend was also 
apparent in Figure 3. As previously discussed, the sample 
processing region are not greatly influenced by the quench block 
temperatures studied. Therefore, several different scenarios for 
cooling the Arc Furnace are possible, including forced air, 
water, or change of phase cooling approaches. As an example, one 
might choose to design a quench block and holder combination 
which would use the conversion of water to steam to remove heat 
and provide a constant reference temperature at the cooling 

block . 

Radial profiles for these numerical calculations are shown 
in Figure 6 a, b, c, and d for the 50, 100, 150, and 200°C, 
respectively, and Figure 7 for the 10 and 38°C simulations. 

Again, these figures remain virtually unchanged in the region 
discussed above. In these figures, the zero line represents the 
centerline temperature, with the solid line representing the 
tungsten surface. The dashed line is the difference between the 
centerline profile and the profile at a radial distance of 2 mm, 
while the additional line is for the difference in the radial 
profile at 1 mm relative to the centerline. The largest 


- 17 - 



RADIAL DISTANCE 



.00 5.0 10 . 15 . 20 . 25 . 30 . 35 . 


AXIAL OISTANCE < mm ) 


Figure 5 . Contour plots obtained for the Arc Furnace with a 

tungsten sample and stainless steel holder. Boundary 
conditions on the external surface of the holder were 
controlled at 10°C (dashed line) and 38°C. (solid line). 
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Fiaure 6. Radial temperature distributions relative to the 
centerline temperature for the Arc Furnace with a 
tungsten sample and stainless steel holder. Boundary 
conditions on the external surface of the holder were 
controlled at a) 50°C, b) 100°C, c) 150°C, and d) 
200°C. 
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Figure 7. Radial temperature distributions relative to the 
centerline temperature for the Arc Furnace with a 
tungsten sample and stainless steel holder. Boundary 
conditions on the external surface of the holder were 
controlled at 10°C (dashed line) and 38°C (solid line) . 
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deviation in the radical profile occurs on either side of the 
sample entrance into the stainless steel holder and quench 
region . 

Figure 8, 9, and 10 and Figures 11, 12, and 13 show a 

similar plot to those discussed above but using water cooled 
aluminum and copper quench blocks/holders, respectively. Figures 
8 and 11 show the axial center line profiles developed with the 
aluminum and copper holders. The maximum axial gradient 
calculated for these materials as holders was 349°C/mm for the 
aluminum holder and 358°C/mm for the copper holders. From an 
examination of the isotherms in Figures 9 and 12, both quench 
blocks will effectively cool the tungsten sample without melting 
the holder material. If one considers the influence of the 
cooling media in Figures 9 and 12, it becomes obvious that higher 
operating temperatures on the cooling media greatly improves the 
uniformity of the thermal gradients in the core of the sample. 
Figures 10 and 13 permit an analysis of the radial temperature 
profiles relative to the centerline temperature. This analysis 
shows that the aluminum holder will provide more uniform quench 
for directional solidification operations. 

3.5 Objective 5. Arc Furnace Performance Verifications. 

The fifth and final objective for research concerned the 
initial performance verification of the S. T. International 
Orbital Arc Welder. The arc welder power supply and 3" orbital 
head were purchased in September 1988. The 3" head was flight 
tested in November 1988. During this testing the "Home switch 
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centerline temperature for the Arc Furnace wicn a 
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Contour plots obtained for the Arc Furnace with a 
tungsten sample and copper holder. Boundary conditions 
on the external surface of the holder were controlled 
at a) 50°C, b) 100°C, c) 150°C, and d) 200 C. 
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Figure c ; nterline temperature for the Arc Furnace with a 
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in the head was found to be defective. The head was returned to 
the manufacturer for repair in December 1988. The manufacturer 
returned the 3" weld head in January 1989 where it was mated with 
the power supply for testing in the laboratory. These tests 
indicated that the power supply was not functional and it was 
also returned to the manufacturer for repair. The repaired unit 
was received in February 1989 and tests were performed to verify 
its performance. The unit passed all tests, welding 3" diameter 

stainless steel and aluminum tubing. 

During the flight test of the 3" weld head, it was noticed 

that the temperature of the weld head body greatly exceeded past 
measurements with 1" weld head. Since the 3" head was ordered 
with a water-cooled feature, it was decided that a system was 
needed to control the head temperature during flight. The system 
had to be light, portable, and dependable. Such a system was 
designed as follows. The basic components needed to cool the head 
included; two Amtrol holding tanks, a Love Model 690-1 
controller, a Hoke 115 volt AC solenoid valve, and associated 
fittings and tubing. The Amtrol holding tanks were selected 
because they contain bladders which can be pressurized for 
operation in microgravity. The cooling water is held in the 
upstream supply tank which is fitted with the solenoid valve and 
pressurized to 75 psia. Hot water is collected in a second tank 
which has its valve stem removed so that its bladder operates 
against a constant back pressure (15 psia) . In operation the 
controller measures the weld head temperature using a type K 
thermocouple and compares it to the set point temperature. If the 
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temperature is higher than the set point, the solenoid valve is 
opened allowing pressurized water to pass through the head and 
into the collection tank. The unit was tested in February 1989 

and performed as expected. 

Additional research was directed at estimating the 
translation rate of the S. T. International Arc Furnace, the 
design of refractory lined collet for improved heat transfer, and 
development of a cooling system to handle the additional heating 
loads generated by the lined collet. Work on the ST 
International Arc furnace included the modification of an 
existing translation mechanism that allowed directional 
melting/solidification of two different samples. Samples of 
aluminum rod and thin walled stainless steel tubing were melted 
and/or heated while translating at various rates to investigate 
the translation speed and power required that would give 
consistent melting zones within the samples. In addition a 
keyboard, color monitor, and color video card were installed in 
the welder power supply to permit easy access to the weld 
parameters . 

To evaluate the potential for the S. T. International 
weldhead to function in a directional solidification mode we 
performed the following experiments. The weldhead was installed 
on an existing jackscrew driven platform designed to move 
horizontally. The weldhead was translated using a variable speed 
DC drive motor connected to the fine thread jacking screw. 

Samples of either 1/2" O.D. aluminum rods or 1/2" O.D. with 1/16" 
wall stainless steel tubing were supported on both ends using an 
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aluminum support stand. The stand was connected to the positive 
terminal of the weldhead. Two inch O.D. aluminum collet lined 
with lava stone of 5/8" I.D were manufactured and served to both 
insulate the sample and contain the argon atmosphere needed 
during arcing operations. An existing color video card was 
exchanged for the monochrome video card used in the IBM-XT which 
is inside of the S.T. International power supply. A Zenith 
keyboard and color monitor were connected to the unit. Addition 
of these components, permitted more ready access of welding 

routines provided with the unit. 

Preliminary tests were performed on stainless steel tubes a^ 
current levels ranging from 50 to 90 amps and translation rates 
of 0.1 to about 5 inch/min. At the highest translation rates and 
maximum power, the sample did not melt uniformly with the melted 
zone resembling a spiral. The best results were obtained for 

translation rates at or below 1 inch/min. 

Tests were also performed on 1/2" O.D. solid aluminum rods 
at conditions similar to those previously discussed. The 
aluminum rod, having a greater thermal conductivity and heat 
capacity, required maximum currents and much slower translation 
rates than did the stainless steel samples. At a translation 
rate of 0.25 inch/min and a current of 90 amps the weld pool only 
extended half way into the sample. Since the sample was not 
preheated, this result suggests that aluminum and, perhaps, large 
rod samples in general, will require preheating prior to the arc 

melting . 
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4.0 CONCLUSIONS AND RECOMMENDATIONS 

The major results of this research show that the current 
configuration of the ARTCOR Model 460-15 Bench Top Tubular 
Laboratory Research Furnace lack sufficient thermal controls to 
provide scientific temperature data and lacks dimensional 
stability to retail the zirconia element in a fixed location. 

Lack of element support on one end is an engineering problems 
that can be addressed with simple modifications of there design 
as shown by using the research version prepared during this 
research . 

The controls problem is a larger problem that requires 
further research. The upgrade of the controller or the inclusion 
of a computer/controller combination that would permit feed-back 
control based on temperature measurements from either 
thermocouples or an optic pyrometer, and would greatly improve 
the research potential of the ARTCOR Model 460-15 zirconia 
furnace research for use with directional solidification in both 
unit gravity and eventually microgravity environment. 

Control techniques should include a direct temperature 
feed-back technique. Analysis or calibration using a constant 
emmisivity optical pyrometer may not predict the processing 
temperature to within limitations placed on directional 
solidification standards. This conclusion is based on the 
spectral radiation models developed for a similar geometry. A 
possible control technique that uses an optical pyrometer 
incorporating spectral radiation within a numerical model to 
adjust the applied power may be possible based on the numerical 
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results. It is apparent from this work, that solving the 
temperature measurement and controls problem is the key to 
advancing this furnace toward directional solidification 
operations . 

Element failure was evidenced with the two elements 
purchased from the manufacturer. Though the manufacturer claims 
that the elements have long life our research did not support 
this claim. This is not to say that all zirconia elements will 
last a short time, it is possible that we received two elements 
from a bad batch, since an experimental element prepared by us 
without input from ARTCOR has operated successfully for over 37 
hour without failure, even after ramped through a thermal 
gradient of over 45°C/min from 900-1600°C. Research with this 
element shows that the basic concepts of using a zirconia ceramic 
tube as a heating source are sound and deserving of further study 
for directional solidification operation under an air atmosphere. 

Both furnace technologies show great potential to advance 
directional solidification to higher operating temperature. The 
Arc Furnace concept should be developed. The Zirconia Furnace 
needs to be further researched to correct the problem found 
during this research, before development. 
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ABSTRACT 


This study focuses on experimental and numerical considerations 
that are required in the design of high temperatures materials 
processing furnaces using zirconia heating elements. 

The experimental study was performed using a commercial 
zirconia tubular research furnace. This furnace was evaluated and its 
application in microgravity experiments assessed. Based on this 
furnace a number of design improvements and modifications were 
suggested. A new zirconia heating element was designed, based on 
oxygen sensor technology and the results from numerical studies that 
were performed in this study. Improvements to the furnace were made 
and together with the newly designed elements received further study. 

One and two dimensional conduction/radiation models for an 
alumina crucible in a zirconia heater/muffle tube enclosing a liquid 
iron sample were solved numerically. Variations in the crucible wall 
thickness were numerically examined. The results showed that the 
temperature profiles within the liquid iron sample were significantly 

affected by the crucible wall thicknesses. 

Experimental and numerical studies were accomplished in 
parallel so that integration of these studies could produce a 
workable model for the furnace while extending the capabilities and 
understanding of the zirconia heating elements. 
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CHAPTER 1 


INTRODUCTION 


This study focused on the development and design of a high 
temperature zirconia furnace. A zirconia tube was used as the main 
heating element. Normally zirconia is an insulator, but when heated to 
temperatures above 1000K it becomes conductive and thus can be used as a 
heating element in such furnaces. The thesis included experimental and 
numerical studies that were performed to understand the heat transfer 
interactions of the furnace. The study is limited to high temperature 
furnaces that used a zirconia ceramic heating element. 

The furnace development was influenced by work presented in the 
literature and the numerical work performed in this study. Numerical 
modeling of a crucible wall enclosing a liquid iron sample with 
radiation/conduction were performed to study the various modes of heat 
transfer. The numerical work was performed to study parameters that 
could affect furnace performance. Among the factors needed to be studied 
in the design of an efficient furnace system were crucible wall 
thicknesses and the transparency of the ceramics used in the 
construction of the furnace. 

In the design of high temperature materials processing furnaces, 
temperatures in the range of 1873K to 2273K (1600°C - 2000°C) are 
anticipated in the hot zones. At these high temperatures, ceramics used 
in the construction of critical furnace components, i.e. muffle tubes 
and crucibles, enter regions in which spectral radiation effects cannot 
be ignored if adequate control of the furnace is to be achieved by open 
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loop control techniques. 

Models used in the open loop control must predict temperatures in 
che samples based on measurements made either at the furnace wall or 
■j/terior crucible surfaces, using either direct measurements or optical 
-vrometry. When spectral radiation properties as functions of wave- 
lengths are considered, the wall assumes a diffuse boundary in which 
some radiative frequencies see an opaque surface at some frequencies and 
a semitransparent surface at other frequencies. 

The scope of this thesis included a review of previous literature 
cn furnace design and development, techniques that are available for the 
measurements of temperatures in semitransparent media, and the numerical 
techniques used in solving combined modes of heat transfer. To 
facilitate the understanding of combined modes of heat transfer and the 
numerical modeling of such processes, a review of the literature is 
given is section 1.1. 

The experimental investigations that were performed in previous 
/ears and a review of the methods used to measure temperatures in 
semitransparent materials is given in section 1.3. The transparency of 
ceramics at high temperatures was a major consideration in the design of 
this furnace, therefore the absorption spectra of solids in general will 
be discussed in section 1.2. 

1.1 ANALYSIS OF COMBINED MODES OF HEAT TRANSFER 

This thesis focuses on how numerical studies could be used in 
conjunction with experimental work to calculate the temperature profiles 
chat would closely resemble real experimental conditions. The 
formulation of the radiation/conduction models in this thesis will take 
into account the radiative properties as functions of wavelengths and 
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temperature . This should closely resemble the actual heat transfer 
processes within the alumina crucible and sample. 

Energy transfer in solids is normally treated as if conduction 
were the only means of transport. This is an erroneous assumption when 
the temperatures are such that thermal radiation is significant and >_he 
material is partially or completely transparent to this radiation. 
Neglecting spectral radiative heat transfer can lead to incorrect 
conclusions and design failure. 

Presently, our understanding of coupled conduction- radiation is 
based upon analysis with emphasis on gases rather than solids (Edwards 
and Matovosian, 1984). Unlike gases, where the conductive flux is small 
and may be justifiably neglected, conduction is always present in a 
solid. Additional phenomena which are not present in gaseous radiative 
transfer must also be considered. According to Kirchoff's laws, the 
emission of radiation is proportional to the square of the refractive 
index. Hence, emission plays a greater role in high index solids than 
gases whose index is approximately unity. There are significant 
differences between gaseous and solid spectral absorption coefficients. 
In contrast to gases whose absorption spectra are established by quantum 
and statistical mechanics, some of the absorption mechanisms of the 
solid state remain to be identified and one must of necessity turn to 

experimental data (Armaly and El-Baz, 1977). 

Little experimental or analytical work directly applicable to 
solids at these extreme conditions has been reported. Difficulties with 
controlling the experimental conditions and accurate measurement of 
temperature at the required levels have curtailed experimental studies. 
For simplicity, analysis is commonly based on unrealistic property 



models and typically concentrated on the radiation dominant situation. 
Simplified procedures for real materials such as those developed for the 
gray model (Olfe, 1963) are required so that the significance of the 

radiative transfer may be assessed. 

There have been many previous studies in combined modes of heat 

transfer. Viskanta and Grosh (1962) solved a one-dimensional 
simultaneous conduction and radiation problem by reducing a non-linear 
integro-differential equation to a nonlinear integral equation. They 
then solved the integral equation iteratively and compared their results 
with Rosseland and Milne - Eddington approximate solutions. Anderson and 
others (1973) further investigated the heat transfer through a plane 
layer of a semitransparent solid experimentally and compared the results 
with the approximate techniques. Using a method similar to Viskanta and 
Grosh, Eryou and Gltcksman (1958) compared the results with their 
experimentally measured one -dimensional temperature profiles and heat 
fluxes within a slab of molten glass. They reported agreements of 
measured values of within 5 K throughout the slab. Grief (1964) studied 
a one -dimensional energy transfer problem of an emitting, absorbing, and 
conducting gas. He considered the thermal conductivity and absorption 
coefficient as having power law dependence on temperature. 

The calculation of radiative transfer in a cylindrical 
nonisothermal emitting- absorbing medium has attracted considerable 
attention in recent years, because of its important applications in many 
high temperature phenomena associated with planetary reentry, nuclear 
explosion, and industrial furnace designs (Kellet 1952. Kingrey . 1976 ) . 
The emerging field of materials processing and crystal growth will pose 
challenging problems in the area of combined heat transfer in 


semitransparent materials. Various 


limits in radiation heat transfer 


problems have been investigated, depending on whether the dominant mode 
of transfer was conduction or radiation (Wang and Tien 1967). 

Modest and Azad (1981) evaluated the radiative heat flux in an 
absorbing, emitting and 1 inear - anisotropically scattering cylindrical 
media. Their study was limited to 1 inear - anisotropy for two reasons. 

First, exact analytical solutions (Lick, 1965) are practically 
impossible to obtain for higher-order scattering phase functions. This 
is due to the fact that the azimuthal direction cannot be eliminated by 
integration as in the planar case, resulting in higher order integral 
equations in the radial and local direction coordinates. Secondly, the 
accuracy of a linear-anisotropic approximation to realistic scattering 
phase functions in the heat transfer applications without collimated 

irradiation had been already established. 

The difficulty in exact calculation of multidimensional radiative 
transfer has prompted researchers to develop approximate schemes such as 
the mean-beam length method (Siegel and Howell, 1980), the zoning 
method, the differential approximation, and the Monte Carlo method 

(Steward and Cannon, 1971). 

Chang and Smith (1970) analyzed the heat transfer in a conducting, 
emitting and absorbing medium bounded by two infinite coaxial 
cylindrical surfaces, for both transient and steady states. They 
formulated the problem in two differential equations (one for the 
radiation potential and the other for the temperature) using 
quasi-steady simplification and Eddington's first approximation for 
radiative transfer. They found that for highly emissive surfaces, the 
interaction of radiation with conduction has a negligible effect on the 
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total heat flux. Ho and Ozisik (1988) used the Galerkin method to obtain 
an exact solution to the radiation part of the problem, and a finite 
difference scheme to solve the conduction part. They solved a steady 
state conduction and radiation heat transfer for an absorbing, emitting 
and isotropically scattering two dimensional rectangular enclosure. They 
examined the effects of the conduction- to- radiation parameter, the 
single scattering albedo, and the aspect ratio on the temperature 
distribution. They considered the boundary surfaces to be black and 
subjected to a constant temperature on one surface and zero temperature 
on others . 

A modified finite difference approach to the problem of transient 
combined conduction and radiation in an absorbing and emitting annular 
medium was presented by Gordaninejad and Francis (1984). They studied 
two different boundary conditions; the first one was a step change in 
temperature on both inner and outer surfaces, while the second type of 
boundary condition was a step flux in the inner surface. The solution 
was shown to closely represent the exact solution for the case of pure 
conduction. Fernandes and Francis (1982) solved a transient combined 
conduction and radiation in a gray absorbing, emitting, and scattering 
medium of cylindrical geometry, using the Galerkin finite element 
method. The results were compared with exact solutions for the extreme 
cases, i.e., pure conduction and radiative equilibrium. 

Some experimental investigations and methods that can be used in 
such studies are reviewed in the section 1.3. 

In the following section, a detailed description of the absorption 
spectra for solids, including some ceramics, will be discussed. 
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1.2 ABSORPTION SPECTRA OF SOLIDS 

In order to understand the semitransparent nature of solids, the 
nature of the absorption spectra of solids with respect to spectral 
radiation at high temperatures needs to be analyzed. 

Radiative transfer analysis requires the spectral absorption 
coefficient which may be obtained from either theoretical models or 
experimental data. For gases the absorption spectra are well known and 
have been reduced to workable models for energy transfer processes 
(Edwards, 1984). Such is not the case for the liquid and solid phases 
where the energy levels are largely determined by the intermolecular 
binding and proximity of molecules. Since these two factors are 
determined by the composition and molecular arrangement as well as a 
number of other factors, the absorption spectrum can no longer be 
obtained from the structure of an isolated molecule and its collisions 
with others. Rather, the entire molecular structure, which may take a 
variety of forms, must be considered to accurately predict the spectral 

absorption coefficient of a solid or liquid. 

The optical properties of solids are thoroughly discussed in the 
literature. All of these are devoted to the interpretation of the 
optical properties from the crystallographic structure and 
spectrographic analysis. Since this area is properly left to the 
expertise of the solid-state physicist the intent here is not to delve 
into the details of absorption but rather to examine some representative 
data for possible simplifications and generalizations which may be used 
for energy transfer purposes. Representative spectral data for the range 
in which thermal radiation is important at temperatures normally 
encountered are presented in Figure 1.1. 
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Common to most materials is a strongly absorbing, broad lattice 
vibration band as represented by the inserted data of Gottlieb (1960) . 

Over the majority of this band the absorption coefficient is quite large 
(10 3 to 10 5 per cm) and the solid is effectively opaque. In the visible 
and near infrared ranges these materials become transparent eliminating 
any interaction between the material and radiation except for reflection 
and refraction. For the purposes of thermal radiation transfer, only the 
semitransparent absorption edge influences the temperature of the 
medium. As with any generalizations there are exceptions such as type II 
diamond which has a narrow weak band that must be considered In its 
entirety . 

Many of the semi transparent regions might be treated empirically 
by an expression of the form - a ^ ^ which is the same as that for a 
gaseous ionization edge (Nelson and Crosbie, 1970). Lithium fluoride Is 
a good example of this empirical fit. Other materials have more complex 
spectra as, for example, synthetic sapphire and silicon. These anomalies 
can be attributed to several resonant modes. Since no general 
correlation can be found which fits any but a limited class of 
materials, experimental data must be correlated separately. 

The spectral absorption coefficient also has a dependence on 
temperature. For lattice vibrations the temperature effects are only 
apparent in the band center (Gottlieb, 1960). To avoid fracturing a 
solid the temperature gradients must be minimal, thereby restricting the 
temperature range normally encountered in practice. 

It can be concluded that for computational purposes a piecewise 
approximation appears to be appropriate for correlating complex spectra. 
Workable comprehensive theoretical models for nonlinear spectral 





behaviors do not appear to be available for all but a few semi- 
transparent solids. 

In conjunction with a heat transfer numerical model, an optical 
pyrometer can be used to calculate a temperature profile that would 
closely resemble real experimental conditions in a furnace. The 
temperature profile can then be used in the feedback control of the 
furnace temperature. High temperatures materials processing furnaces 
that are used for zone refining and other processes, depend on accurate 
control of their temperatures for successful processes. 

The formulation of the numerical model which uses the spectral 
properties for the materials is discussed in chapter 2. 

1.3 EXPERIMENTAL INVESTIGATIONS AND METHODS 

A long range goal of this research is to show that it is possible 
to deduce an accurate temperature profile from measurements made by an 
optical pyrometer if these measurements could be used in conjunction 
with a heat transfer model of the crucible. 

Many applications are encountered in applied physics and 
engineering where it is necessary to measure the temperature 
distribution within a semitransparent material at high temperature. 
Examples of such materials are combustion products in a combustion 
chamber, glass when solid or molten, crystals, plasmas, etc. It is 
difficult to measure the temperature distribution of such materials. 
Probes cannot be readily used because they disturb the temperature 
measurements and radiation fields, and thus produce erroneous results. 
Also probes must be specially designed for high temperature applications 
to keep them from melting; for example pulsating thermocouples for 
combustion applications. An alternative approach is to remotely measure 
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che radiation field emerging from the medium and recover the temperature 
distribution from these data (Chupp and Viskanta, 1974). This is 
possible because the emerging radiation field from a semitransparent 
medium is a function of the temperature distribution. For the limiting 
case of opaque materials, the remote sensing approach reduces to 
conventional pyrometry and only a surface temperature is recovered. 

Thermal remote sensing specifically involves using a spectral 
sensor (spectrometer or monochrometer with detector, radiometer with 
filters, etc.). This equipment is expensive and may not be required if 
accurate modeling of the semitransparent materials can be achieved. 

A number of experiments have been conducted for gaseous, liquid 
and solid materials with radiative transfer. Hany of those are of an 
applied nature undertaken for purposes other than identifying the role 
of radiative transfer and its coupling with other energy transfer modes. 
Since the experimental methods developed in these works are applicable 
co a more fundamental objective, they are included in this study. Most 
of the techniques which will be discussed are applicable to all 
materials and are independent of the material phase. Therefore this 
review will not concentrate on the phase of the material. The different 
types of techniques, their limitations, and applications in temperature 
measurements will be discussed in the following sections. 

1.3.1 Thermometry Methods 

Several investigators have applied conventional thermocouple end 
resistance thermometer probes for the measurement of temperature ulth 
significant radiative transfer. The presence of the probe alters 
temperature field by partial blockage of the radiation, conduction along 
Che probe, probe mount and leads, and the absorption and reemis 
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radiant energy by the probe. Since all of these errors depend on the 
physical size of the probe, they can be minimized by appropriate design. 

Nishimura and others (1968), examined the validity of a gray 
analysis for simultaneous energy transfer in glass. It was found that 
this model failed to predict the experimental results. It is of interest 
to note that the gray assumption overpredicts the temperature for a 
solid whereas it underpredicts that of a gas. This behavior can be 
attributed to the fact that solids have but a few regions of high 
spectral bandpasses whereas gases are transparent over a large portion 
of the spectrum. Fowle and others (1969). experimentally studied the 
unsteady heat transfer through glass and measured the surface 
temperatures with thin- film gold resistance thermometers. 

The effect of subsurface temperature gradients on the total 
emittance of a translucent solid was examined by Cox (1963) using small 
blackbody cavities drilled into the specimen normal to the gradient. 
Internal temperature measurements were made by focusing an optical 
pyrometer on the aperture of the cavities. The validity of this approach 
is questionable since the cavities appreciably disturb the radiation and 

conduction fields. 

1.3.2 Optical Methods 

Optical interference methods which have been widely used in heat 
transfer studies eliminate several errors since no external probes are 
introduced into the material. Also, they have the advantage that the 
entire field may be observed instantaneously. The earliest application 
of an interferometric technique to a situation in which radiation is 
significant appears to be that of Gill and Goody (1964) , who measured 
the critical Rayleigh number of a model atmosphere heated from below. 
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* Michelson Interferometer which passes the beam through the test madlum 
twice was used to obtain the temperature distribution. The investigators 
were very meticulous in the design of the experiment and error analysis 
although the correlation with analysis is poor because of the 
statistical nature of the convective instabilities. Ho account of the 
refractive error which the Michelson interferometer accentuates was 

given. 

For solids the temperature field may be Indirectly measured with 
photoelastic optical methods (Veduta and Leontlvlch. 1967), which 
measure the thermal stress. This approach has had little application 
because of the limited accuracy inherent in correlating the 

birefringence patterns and temperature. 

The variation with temperature in the transmittance of molten 
glass at 0.6328 pm was employed by Eryou and Cllcksman (1972). This 
method is more practical than an interferometer in that it is simpler to 
us, and less sensitive to external error sources such as vibration and 
convection. Since the transmittance as a function of temperature is 
many cases unknown, elaborate calibration techniques were required. 

1.3.3 Inversion Methods 

The inference of the Internal temperature from measurements of 
radiant energy emerging fro. the surface (Inversion) has the advantage 
that no probe, or external stimulus need be applied to disturb the 
fields The earliest attempt using inversion technique, to measure the 
temperature was made by Beattie and Coen (1960). They recognised that 
the spectral radiation emitted from a semitransparent medium is 
dependent on the Internal as well a, surface emission. A pyrometer 
capable of remotely measuring the emitted radiation in two wavelength 



bands, a semitransparent band to measure a weighted internal temperature 
and an opaque band to measure the surface temperature, were developed, 

The logical extension is to increase the number of spectral bands 
in order to resolve the entire temperature profile. This extended 
approach has been of considerable interest for astrophysical and 
meteorological applications. For the atmosphere which can be modeled as 
a plane layer, the emergent radiation is a weighted integration of the 
local composition and temperature. Inferring the temperature then 
requires the inversion of this integral equation. 

In chapter 3 experimental investigations that were performed using 
a commercially available AUTCOR Model 460-15 Bench Top Tubular 
Laboratory Research Furnace are discussed. The study also considers 
combined mode spectral radiation- conduction modeling of semitransparent 
materials using finite differencing techniques. 

In chapter 4, the results are discussed, with conclusions given In 
chapter 5. 


CHAPTER 2 


FORMULATION AND ANALYSIS OF SPECTRAL 
RADIATION EFFECTS IN A ZIRCONIA FURNACE 

The understanding of the heat transfer processes within an 
alumina crucible containing a liquid Iron sample Is an Integral part 
in the design of this furnace. Heat transfer modeling at these very 
high temperatures is complicated by the semitransparent nature of the 
materials. The Interaction between conduction and radiation 
the crucible wall also needs to be understood for an efficient design 

of the crucible wall thicknesses. 

in this chapter a general steady state analysis of coupled 
conduct ion/rad la cion heat transfer analysis through semi-transparent 
media will be developed. In section 2.1 a on.-dl«.nslon.l spectral 
r.diaclon/conductlon model In a semi-transparent material was 
developed for a zlrconla furnace, operating at 1873K (1600 C) , with 
alumina crucible containing molten Iron. This model was developed 
to understand how the variations In the radiative functions affected 

the temperature profiles within the system. 

The one dimensional radiation conduction model was able to 
predict the temperature profiles within the alumina crucible wall. 

The model was also used to verify the solution technique, that were 
used. The model could not predict how the variations In wall 
thicknesses influenced the temperature profiles within the crucible 
wall. Therefore in order to solve this problem a two dimensional 
model had to be developed. In section 2.2 a two dimensional spectral 
model for the semitransparent material was developed. 
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Both models included variable radiation properties namely, 
transmissivity, emissivity, and absorptivity all as functions of 
temperature and wavelengths. Numerical results for the temperature 
profiles were obtained using finite difference techniques. 

2.1 ONE DIMENSIONAL ANALYSIS 

A one dimensional spectral radiation heat transfer model for a 
semitransparent alumina crucible containing liquid iron sample was 
developed. The radiation/conduction model included emissivity, 
transmissivity, reflectivity, and absorptivity that are both 
functions of wavelengths and temperatures, as discussed in section 

2 . 1 . 2 . 

2.1.1 One Dimensional Physical Model 

Figure 2.1 shows the cross section of a cylindrical alumina 

crucible containing liquid iron and enclosed by a heated furnace 
wall. The furnace wall temperature is controlled so that the outside 
surface of the alumina crucible wall (a) is maintained at 1873K 
(1600°C). The problem is modeled in one dimension as shown in Figure 
2.2. In Figure 2.2, region II, the space between the muffle tube and 
the alumina crucible, was assumed to be at vacuum. The alumina 
crucible is considered to be an absorbing, emitting, and isotropic 
scattering media, with both conduction and radiation taking place. 

The alumina crucible will be considered as semitransparent in some 
regions of wavelengths as dictated by its transmissivity (case 1) and 
compared to a totally opaque wall (case 2) which will be used as a 
reference condition. Thermal conductivity will vary with temperature 
and the modeling reflects this. Radiation properties will be 
considered to vary with both temperature and wavelength. Finally, in 
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Furnace Wall 



Figure 2.1 A Cross Seccion of a Cylindrical Alumina Crucible 

Containing Liquid Iron and Enclosed by a Heated Furnace 
Wall. 
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region IV, we .sinned liquid Iron with conduction In the radial 
direction. Convection effects were Ignored In the liquid Iron sample 
as It is assumed the unit will operate In microgravity . Finite 
difference techniques were employed to solve this problem using the 
coordinate system shown in Figure 2 . 2 and discussed in section 2.1.3. 

2.1.2 Physical Properties 

This section details the various radiative and physical 
properties of alumina and liquid Iron sample. The spectral radiative 
property functions were obtained by correlating materials property 
data using high order polynomials. The objective was to develop a 
combined radiation and conduction heat transfer model for an alumina 
crucible in thermal cohtact with a molten iron sample In order to 
demonstrate the importance of spectral modeling. 

Figure 2.3 shows the e.lsslvity of alumina a. • function of 
temperature, while Figure 2.4 shows the emls.ivlty a. a function of 
wavelength at 1873K. Figure 2.3 shows the transmissivity of alumina 
as a function of wavelength over the rang, of 1873 to 1973*. Not. 
that at 1873* approximately 80% of the transmitted energy for alumina 
occurs between 2 to 6 microns. It is for this reason that 
transmissivity cannot be ignored in high temperature radiation 
problems. The thermal conductivity of alumina as shown In Figure 2.6 
changes by almost 7e fro. 1773 to 1873*. Finally. Figure 2.7 shows 
the thermal conductivity for liquid iron, which changes by about 30» 
of its value at 1773* to its value at 2000*. This represents a 
significant error if it were to be assumed constant over the 
temperature range examined here, Toulohian (1967). 

This study considers all the radiative properties, namely 
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Figure 2.7 Thermal Conductivity of Liquid Iron as a Function of 
Temperature . 



transmissivity, absorptivity, and a.lsslvlty as function, of both 
■wavelengths and temperatures. This Is Important since at these 
such high temperatures (1773 - 2273K) the radiation properties vary 
due the transparency of the ceramic. Therefor, the transparency of 
ceramics at these high temperatures is important in model development 
and must be included if ve are to advance high cemperacure furnace 

technology . 

2.1.3 Governing Equations and Boundary Conditions 

SLnce only conduction is considered in region IV (see Figure 
2.2), the thermal gradient in region IV will only be a function of 
the heat flux established in region III. Thus, the conduction in 
region IV will be determined from the heat flux at the Interface (b) 
after the results in region III are converged to a solution. That is, 
we shall solve the conduction and radiation within the alumina 
crucible by considering the interface with the liquid iron as a 
boundary condition with constant heat flux at steady state. 

The boundary conditions at the surface of the alumina crucible 
for the two cases studied were as follows. For the semitransparent 
crucible (case 1) we have emission, absorbed irradiation, and 
molecular energy transfer between the interface and the heated 
furnace wall, thus the energy balance is as given in equation (2.1): 

-k ^ - h 1 (T 1 - Tq ] + < 0 pq (e bAopq (T l ) * •bXopq^O^ (2 ' U 

dx x-0 

For the opaque reference condition (case 2), -he temperature 
surface of the crucible is considered constant at 1873 K. 

At steady state in the absence of internal heat sources, 
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conservation of energy requires that 


( 2 . 2 ) 


d dT 

_ -k(T) — + Qr ( x) 

dx 

• 1 • « that the sum of the conductive and 

»• — £Ud ' 11,15 h at nux The radiative flu* 

radiative flux must equal the constant 

ov .r the entire rang, of wavelengths Is therefor. 

V „ - U 

o 0 

s , nrensiev of radiation in 
T + (x) an d I*>(x) are the spectral 

where I an * _hvp1v As 

, . (ll < directions, respectively, 

the forward (w > 0) and backward U < 

of radiation, the intensity is coup 
a result of Internal emission of radlatl 

to the temperature distribution in the solid. Now insider te 
equations of transfer and following the work of Sparrow and 

(1970) , we have 


dI + A 

«bA^ 

- I + A 1 

(2.4a) 

dx 





< bi (x) 


(2.4b) 

dx 

* * 




for 0 < w < I- Hot. that I*g corresponds to positive values 
I* is associated with negative values of 
LI time being let the boundary conditions for equations 

(2.4a) and (2.4b) be 


I + A (r x . 4) 


0 


(2.5a) 


- I+ x (0, M) 


r x (r x , M) - 1 X (r 0 V ^ 


r X ■ r 01 


(2.5b) 


by defining 

/** , ^ (2.6) 

r x - j *^ x > dx 

0 

called the optical thickness. By using either an integrating factor 
or the method of variation of parameters, the solutions to equations 
(2.4a) and (2.4b) are found to be 

l\ (r x . M> - 1*1 <0>.-V" 

a 1/, f *<. ;k « b »(t>l«' (r »' c)/ ’‘ 

0 

0 < H < l (2.7a) 


I‘l (e t . „) - I‘l (toi)« (r0i '' i>/M 

j 1*1 


- 1/* 


e bx U) ]e X 


-(r c)/m 


dt/M 


tx 


-1 < M < 0 (2.7b) 


Equation (2.7) describes the radiation field in terms of the 
temperature field within the medium through Planck's constant e bx 
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Subscicuclng equations (2.7a) and <2. 7b) into equation (2.3). 
:b, expression for rhe on.-dinensional coral radt.cion flux becomes 


Q R (r A ) 2 *fl\ (O' " d " 

0 o 

•/r» (-ov-'"''' 01 '' 117 ' ' iB 


2 y t «i e bA (c)E i (r ^ ■ t)dc 
0 

■ 2 J { «*« b j,U>E 2 (t * r ^ )dC ]dX 


( 2 . 8 ) 


where 


VO d " - / 


dy 


(2.9) 


are the exponential Integral functions. Differentiating equation 
(2.8) we gee 




' 2 * [ 


d Q R (r x )/dr x -/l «»**«b^ I V 

0 


/-■ 


(0, M) 


e' r X /M dM 


* p - X (r 0i .-o)*‘ C, 0 1 " ,i>/, ‘ 

o 

• 2_/"<X*bl(0 E l ( f 1 ' >,<U <2 ' l0> 


Equation (2.10) therefor, describe, the divergence of the radiation 

cerm in the energy equation (2.2). 

. / t in') into the conservation of energy 
Substituting equation (Z.LU) into 


equation (2.2) we get 
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v - m/.\ to. - 


/' 


l\ (r 0 A ,-M)e" (r OA' r > )/M d M ] 


f * k 

■ 2 [ ^ A e bA (t)E 1 ( | r A - c [ )dt ) dA (2.11) 


Assuming diffuse surfaces, the intensities may be expressed in terms 
of radiosicies by the following 


I + a (0,m) - B 1a /» , and IV'OA'^ “ Z 2\ /n 


( 2 . 12 ) 


Substituting equation (2.12) into equation (2.11) 


-f 


d/dx[k(dT/dx) ] -/( <** x e bA (r) ' 2B lA E 2 (r A> 


+ 2B 2A E 2 (r 0A “ r A ) 


r k 

2/ t *i«bX<'> E l< I’-x • c| )dt 1 dX 


(2.13) 


where the radiosities B nA 's are defined by 


B 1A “ f lA e blA + 2(1 ' ( 1A > ( B 2A E 3 (r OA ) 


+ J ^ X e bA (t)E 2 (t)dt) 


(2.14a) 
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&2A “ '2A e b2X 


+ 2(1 - <2A } { B 1A E 3 (t 01 ) 


/* * 

+ / ' c A e bX < ' C ^ E 2 ( ' r 0A * c)dc ' 


(2. l^b) 


where solving for Z 2 \ we 8 eC 


2(1 - 


>2X 


2E 3 (1 - *i\) + { 2X e b2X 


l *‘U e bU 


- 2(1 


- f 

0 

:r ox ) j * 


' c i e bi (c)E 2 (C)dt 


Using 


l/E 3 (ro x ) / *x e bX (OE 2 (r 0X 


dimensionless variables previously defined 


r n \ * C)dc) (2.15) 


T X - J K x (x)dx 


e - T/T 0 


(2.16) 


equation (2.13) becomes 


m 

c 2 x T 0 d/dr [ k (dS/dr) 1 - A - 2B u E 2 (r x ) 


+ 2B2^E 2 ( f ox “ f X^ 


2 f « i <«e‘‘T 0 4 E 1 (|' J - t|)dt ' dJ 


(2.17) 


In equation <2.17,. the first tar. on the right-hand Sid. represents 
the emission fro. th, -dim. .t fg. the second two terms represent 
tne radiation absorbed at point , 4 . end th. fourth term represents 
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the emission from all the other points in the medium absorbed at 
point r 

2.1.4 Numerical Formulation 

A finite difference technique was used to solve equation 

(2.17) , as [A] [9] - [B] . The coefficient matrix [A] was composed of 
the conduction terras with variable thermal conductivities, while the 
force vector [B] was composed of the right-hand side of equation 

(2.17) multiplied by T 0 V*v Since the thermal conductivity varies 
with temperature, we can write the left-hand side of equation (2.17) 

as 

dq/dr - R.H.S. of eqn (2.17) (2.18) 


where q - -k d9/dr . 


The finite difference formula for the above equation is 


[ k(0 n+ i/2 ) (®n+l * V - k < e n-l/2> < e n * Vl 5 ^ ' R,H,S< 

of eqn (2.17) 


__ , /<a , . ft + Vc i) our coefficient matrix 

By defining k(0 n+1 / 2 ) “ °- 5 < lt n K n+1 ; ' 

consists of 


a n ’ 

n- 1 

- 0.5(k n . 1 


a n’ n 

--(0.5k n .i 

♦ k n - 0. 5k n+ i 

a n- 

n+1 

- 0.5(k n ♦ 

k n-l )/h2 


)/V 


(2 . 20a) 
(2.20b) 
(2.20c) 
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propercl.es, was found to be insufficient to study the effects of 
variations in crucible wall thicknesses. Variations of temperature 
profiles within the crucible wall are required to determine the 
optimum crucible wall thickness for samples containing molten iron. 

To solve this problem, a two-dimensional model was developed to 
improve the overall understanding of heat transfer interactions 
occurring within the zirconia furnace. The two-dimensional analysis 
more closely models both the spectral radiation and conduction 
occurring within the zirconia furnace. This model was used to 
determine the temperature profiles within the crucible wall and the 
liquid iron sample. The knowledge of the temperature profiles is 
important in the design of this furnace, since, as will be shown, the 
location of the maximum in the temperature profile was a function of 
the crucible wall thickness. 

An important requirement in the design of an efficient furnace 
is that the maximum temperature location be at a minimum distance 
from the crucible/liquid iron interface, resulting in higher 
temperatures within the sample. If the maximum temperature within the 
crucible wall occurs near the crucible/liquid iron interface then the 
crucible wall can basically be neglected during control of the 
zirconia furnace. 

In this study, as with the one dimensional model, we have 
considered all the radiative properties, namely transmissivity, 
absorptivity, and emissivity as functions of both wavelengths and 
temperatures . 

2.2.1 Two Dimensional Physical Hodel 

The model considered here is shown in Figure 2.8. A heating 
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An outline of che alumina crucible containing liquid iron 
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Figure 2.8 
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2.2.2 Governing Equation* 

Ac steady scace in Che absence of internal heat sources, 
conservation of energy requires that 

-V . ( 7T ) + 7 . Q r - 0 (2.21) 

be satisfied. This implies that the sum of the conductive and 
radiative fluxes must equal the total heat flux. For an annular 
geometry with cylindrical symmetry, equation (2.21) in cylindrical 
coordinates becomes 

- [ l/r 3/dr (k 3T/3r)]- 3/3z(k 3T/3z)+ 3Q R /3r - 0 (2.22) 

where the first two terms represent conduction in the radial and 
axial directions respectively. The second term represents the 
radiation in the radial direction. To partially simplify the 
analysis, radiation is assumed negligible in the axial direction, 
within the crucible wall. Also in the liquid iron, only conduction is 
considered to be the dominating mode of heat transfer. 

Using dimensionless variables defined as follows 

m 

, -j* x r<U 

e - T/T ref 

!7 - Z/Zq 

equation (2.22) becomes 


(2.23) 

(2.24) 

(2.25) 
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13 ae 1 3 30 I 3Q R (r) 

r dr dr x 2 Z 2 q drj 3q N dr 


(2.26) 


where we define the radiation conduction parameter N as follows 


N - « 2 /^ s T 3 «f 


(2.27) 


The transport equation, using the coordinates shown In Figure 
2.8 for a material with refractive index n, Is given by 


dl ^/d r 


P\ l \ + l bX + <J x G A (r)/4 W 


(2 . 28) 


where 1^, the intensity and is composed of I + and I* which are 
dependent on the direction of radiation, and G^(r) is the integrated 
intensity defined below. 


»/2 2 * 

C A (r) - 2 If I^C0Sad7da 
0 0 


(2 . 29) 


Therefore, the net radiative flux directed radially inward is given 
by 


q A (r) 


*/l 2» 



COS 2 ocCOS 7 d7d« 


o 0 


(2.30) 


From the above equation (2.30), the expression for the radiative flux 
directed inward in nondimens ional form is 


w/2 *i.n‘ 


Qd ( T 


) - Ur^nr j ( r , (3) - *f( r , a 3) ) cos 2 «COsfld0d« 


r vailin' 

* Ur^xr l* 2 *(r.*.0) - ♦ 2 “< f ' a ^ ) ’ “ COS 2 «COS3d3d« 

0 .tn' l (f 0 / ri ) 

where * is the nondimens ional ized radiation potential defined as 


(2.31) 


* - it 1/ n 2 (7 s T 4 ref 


(2.32) 


and I is the intensity of radiation. 

Substituting equation (2.31) into equation (2.26) gives the final 
form of the governing equation. 

2.2.3 Boundary Conditions 

Again, we will consider two kinds of boundary conditions. 

First, case (1), we shall assume that the surface of the crucible is 
opaque. Secondly, case (2). this same surface will be considered as 
semitransparent for some regions of radiation wavelengths. The other 
boundary condition for both cases will be a zero heat flux at the 
center of the liquid iron sample. This assumption is justified by the 

cylindrical symmetry of the model. 

Since the model now has two separate regions where heat 
transfer is being considered, that is region II and III. there exists 
another restriction at the interface of the crucible wall and the 
liquid iron sample. At this interface between the crucible and the 
liquid iron sample, we shall assume that the temperature is 
continuous, thus the heat fluxes on both sides will be identically 
matched. This assumption is justified by the fact that the crucible 

in thermal contact. If we assume 


wall and the liquid iron sample are 
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zero contact resistance, then the temperatures at both sides of the 
interface are equal. 

For case (1), the opaque case, the temperature at the crucible 
external wall will be calculated from the radiating zirconia element. 
This implies that at the boundary r - 1 will be fixed at a known 
temperature. For the semitransparent case (2), the boundary condition 
at the surface of the wall is given below 


-k dQ/dr 


r-R 


( opq t e bAopq ^1 ^ 


e biopq^0^ 1 ® 


(2.33) 


When we consider the interface, as mentioned earlier, the 
temperatures at the crucible side and the temperature on the liquid 
iron side have to be identical. Therefore we can write 


^crucible ” ^iron aC ^ " ^interface ^ 

For both cases the boundary condition at the center of the liquid 
metal sample is that the temperature gradient is zero. Thus we have 
at r - 0 , 


d6/d r - 0 


(2.35) 


Now for the iron sample we have two more boundary conditions for the 
axial (z) direction. Using the coordinate system applied to the model 
(see Figure 2.8), at the center of the sample, ( r) -0), we have 


dS/dr; - 0 


(2.36) 
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while at each end of the crucible wall (r, - l.-l). «• have 


d9/drj - 0 . 5 Q r 


(2.37a) 


de/d rj - -0.5Q R 


(2.37b) 


Equations (2.26) and (2.28) toge 
the final equations needed Co solve the problem 

2.2.4 Numerical Formulation 

A finite difference formula 

, r » i r a 1 - f radiation terms 
is given by [ A 1 l 9 I 1 

[ a ) consists of the following 


ther with the boundary conditions are 


cion of the energy equation (2.26) 

The coefficient matrix 


a n-l.m ‘ k n-l/2/S 


(2.38a) 




a n+l , m “ k n+l/2/8 r m 


2 3 2 

a n , m - 1 *n-l/2 ' u 


V~1 - k n*l/2 /h2 “’ z2 ° 


is defined as 0.5(k n + k n+l^ and k n-l/2 


(2.38c) 


(2 . 38d) 


(2 . 38e) 


is defined as 


where k n+1/2 
0.5(k n . 1 + k n ) with k n being the 

a a - At and h - &V • 

°n' 6 


thermal conductivity at temperatur, 



The left hand side coefficient matrix represents the conduction 
terms (equation 2.26), while the right hand side represents the 
integral radiation terms in equation (2.31). The finite difference 
equations for the conduction part included two spatial variables, 
namely r and n in the radial and axial directions respectively. For 

the crucible wall, the integrals on the right hand side of equation 

(2.26) were calculated in the following manner. For an assumed 
temperatures, the integrals were calculated using simpson's method 
for the range of wavelengths given between each of the nodes. This 
provided the total radiation properties at those nodes. Since this 
problems represents a double integration of complex mathematics, 
simpson’s rule was first used to integrate over discrete wavelengths 
to reduce the problem to a single integral of sums of the integrated 
wavelengths, which were further integrated using a three point 
gaussian quadrature. This calculation provided the elements for the 
force vector at each of the nodes. For this force vector the 

temperature was assumed linear between the nodes. 

An implicit method was used to solve the equations in region II 
(the crucible wall), and the flux and temperature leaving the end 
node of the wall was used to start the calculations in region III 
(the iron sample). When both the solutions converged, a test was 
performed to compare the fluxes at the interface. When die thermal 
fluxes were within the prescribed error of 10' 2 the profiles were 


assumed to converge. 


The results for both the one and two dimensional models are 
discussed in chapter 4. The model results were used in conjunction 
with experimental results to modify the zirconia furnace heating 



u2 


element:. Appendix (B) and (C) show the computer code developed to 
solve this problem. 


CHAPTER 3 


EXPERIMENTS AND PROCEDURES 

This chests addressed improvements In Che design of high 
temperature furnaces that use zirconia ceramic heating elements. 

The goal was to adapt, Improve, and redesign a zirconia furnace for 
use In materials processing experiments. This chapter will discuss 
the experimental apparatuses and procedures used to examine both a 
commercially available zirconia furnace and a furnace developed 
during this research. 

Both furnaces, used a zirconia heating element that can 
operate continuously at temperatures above 1873K (1600°C) and can 
reach limits of 2273K (2000°C) . The above temperature limits were 
based on the mechanical deterioration of the elements and 
insulating ceramics, the atmospheric air operating environment and 
electrical damage resulting from current flow within the zirconia 
heating element as the resistance decreased with increasing 
temperature. As mentioned earlier, the zirconia element is 
basically an insulator at temperatures below 1073K (800°C) but 
becomes a conductor above this temperature. As the temperature is 
increased, the resistance of the zirconia element is reduced as 
will be shown in later sections. 

The major components for each furnace consisted of a main 
furnace assembly and a number of auxiliary systems that are 
necessary to perform the experiments. Section 3.1 describes in 
detail both the commercial and research furnaces and auxiliary 
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subsystems. The auxiliary subsystems included such items as 
temperature measurement devices, control and power supply system, 
and fluid delivery systems. This section also highlights the design 
improvements that were made in the furnace arrangement to include 

an alumina crucible assembly. 

To develop this furnace for eventual application to 
directional materials solidification, several experiments had to be 
designed to verify its performance. Initial experiments were 
performed using an ARTCOR model 460-15 Bench Top Tubular Laboratory 
Research Furnace. This furnace will be referred to as the 
commercial ztrconia furnace (CZF) . The commercial furnace and its 
associated auxiliary subsystems are discussed in sections 3.1.1 and 
3.1.2, respectively. The procedures for this furnace are discussed 
in section 3.2.1. Using the CZF, preliminary experiments were 
accomplished to evaluate its performance. Results from this furnace 
formed a basis on which new designs and procedural modifications 
were made which ultimately led to an improved furnace. This 
improved furnace will be referred to as the research zirconia 
furnace (RZF) . The equipment and procedures used with the research 
furnace are detailed in sections 3.1.4 and 3.2.2. 

A parallel study was performed on the numerical modeling of 
semitransparent materials heat transfer. The study assisted in the 
overall understanding of the heat transfer within the zirconia 
furnace, of the drawbacks of optical pyrometry and the thermal 
profiles within the crucible wall. The information obtained from 
the numerical model was used to select the internal components of 
the furnace that is detailed in section 3.1.3. 
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3.1 DESCRIPTION OF EXPERIMENTS 

This section describes the furnace assembly. A description of 
the ARTCOR Model 460-15 is given in section 3.1.1 and modifications 
co the furnace assembly are then discussed in section 3.1.3. This 
section also describes the various components of the experimental 
setup for the two sets of experiments, set *1 and *2 . The setup for 
the commercial zirconia furnace experiments (set *1) are given in 
section 3.1.2 while the setup for the research zirconia furnace 
experiments (set #2) are given in section 3.1.4. 

A number of auxiliary modules had to be integrated into the 
experiment so that successful investigations could be accomplished. 
An infrared pyrometer (Omega type OS-1000-HT), with a water jacket, 
was included for temperature measurements. Specifications of the 
pyrometer are given in Appendix A. A stand was required to support 
the furnace, optical pyrometer, and ventilator. The stand was 
designed to protect the furnace from vibrations and was fitted with 
clamps to hold accessories needed during the experiments. 

Additional procedures were developed after thorough 
investigation of the commercial furnace operations. 

3.1.1 The Commercial Furnace Assembly 

The furnace assembly is the core of the experimental set up. 
Figure 3.1 shows the commercial furnace arrangement for the ARTCOR 
model 460-15 bench top tubular laboratory research furnace. The 
main heacer element is a zirconia tube. The zirconia acts as a 
conducting material at high temperatures, thus acting as a heater. 
Lead wires are actached to the zirconia tube at the cop and bottom. 
The zirconia heacing element is loosely held in place by cwo 
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heater element support*; an upper support, and a lower support. 

Both supports are made of low density alumina ceramic tubing. At 
the middle of the zirconia element is an insulation assembly, The 
insulation cylindrical blocks consist of high and low temperature 
zirconia insulating blocks. These insulating blocks are held in 
place by highly porous ceramic rings, labeled the Insulation 
support assembly (see Figure 3.1). The above forms the main heating 
element assembly. 

The zirconia heating element and the insulating blocks, 
including all the support rings, are contained in a preheater 
module. The preheater also acts as an insulator after the zirconia 
main heater element has ignited. The assembly also consists of a 
pressed kaolin top and bottom and an outside casing made of 
anodized aluminum. 

Electrical terminals that feed the current to both the 
preheater and the zirconia main heater element exit from the top of 
the assembly. The furnace assembly is attached to its control and 
power box by two protruding pipes. The whole assembly sits on a 
pressed ceramic board. 

3.1.2 Experimental Layout for Commercial Furnace 

The experimental layout consists of the furnace arrangement, 
temperature measurement devices, the fluid delivery systems, and 
the control and power supply. The furnace assembly is the major 
equipment for this experiment and has been detailed above. In order 
to adapt the furnace to material processing experiments, a number 
of auxiliary systems were installed. This section will discuss 
these auxiliary systems. The experimental layout is shown in 


Figure 3.2 This set up was used to perform experiment set -1. 

For the initial experiments two temperature measurement 
devices were installed to permit temperature measurements on the 
commercial unit. A type K thermocouple with a digital display was 
used to monitor the centerline temperature within the heating 
element. An Infrared pyrometer was needed to monitor the 
temperature of the furnace via a viewing hole provided with the 
unit. This pyrometer measured the temperature of the element using 
an average emissive power of the element, therefore the emissivity 
of the zirconia had to be manually adjusted in the pyrometer 
electronics package to reference the pyrometer. A goal of the 
numerical portion of this research addresses the accuracy of using 
an "average- emissivity. A second goal of this research was to 
correlate the pyrometer readings and numerical results with those 
of the centerline temperatures within the alumina crucible. 

The fluid delivery system cooled the pyrometer at high 

operating temperatures. The pyrometer was equipped with a 

_ a 4 fc < nl pt lens The water temperature in 

water-cooled jacket around Its Inlet Lens. 

the Jacket was normally kept around 10 K above room temperature. 

The control and power supply to the main heater and preheater 
formed the fourth subsystem. An Eurotherm controller (Model 818P) 
was used to supply the power to the heating elements by switching 
the channels allocated to the preheater and main heater elements. 
This switching presented some problems in the furnace operation, as 

will be discussed in later sections. 

AS discussed above, the ARTCOR Model 160-15 zirconia furnace 

was equipped with the necessary auxiliary systems so that an 



£ure 3.2 Layout of the Experiment including the Auxiliary Systems 
for Experiment Set *1. 






evaluation on its performance could be made. As modifications were 
made to the commercial furnace assembly, further modifications to 
the auxiliary systems were needed. These changes formed the setup 
for the zirconia research furnace experiments (set #2) . This set up 
is discussed in following sections. 

3.1.3 The Research Furnace Assembly 

The design changes in the furnace assembly are shown in 
Figure 3.3. To adapt the furnace for materials processing 
experiments, alumina crucibles and supporting ceramics were 
installed. These components were used to form the muffle tube and 
support small diameter crucibles which could be filled with pure 
component metals with known melting points. These experiments were 
designed to calibrate the optical pyrometer, evaluate the pyrometer 
in closed loop control schemes, and test the numerical models 
developed for semitransparent alumina crucibles which were 
discussed in chapter 2. 

One part of the commercial furnace was redesigned. The roof 
of the furnace was redesigned from a 1" thick by 6" diameter yttria 
stabilized zirconia blank to partially restrain the element in both 
the radial and axial directions. This restraint would be needed to 
operate the furnace in a microgravity environment. The redesigned 
part also decreased heat losses from the top of the furnace while 
supporting and stabilizing the main heating element. 

The main heating element, was replaced by an experimental 
yttria stabilized zirconia heating element. This new element is 
discussed in section 3.2.4. A type K thermocouple was Installed 
between the main heater and preheater to accommodate a preheater 









power supply arrangement as discussed in section 3.1.2. 

Alumina crucibles and supporting ceramics were used to form a 
muffle tube and support small diameter crucibles which could be 
filled with pure component metals with known melting points. These 
experiments were designed to resemble process conditions that exist 
in materials processing environments. In addition a type S 
thermocouple with electronic readout was added to measure the 
centerline temperature immediately adjacent to the viewing hole. 

With the furnace equipped as shown in Figure 3.3, experiments 
were performed to compare the thermocouple (centerline) temperature 
measurements to the optical pyrometer readings from the external 
wall of the zirconia furnace. 

3.1.4 Experimental Layout for the Research Furnace 

Figure 3.4 shows the experimental layout of the research 
furnace and its auxiliary systems. This setup was used to perform 
the experiment set #2. In addition to the pyrometer, a type S 
thermocouple with a digital display was installed to monitor the 
centerline temperature in the crucible. These measurements were 
needed to correlate the pyrometer readings with the centerline 
temperatures . 

In addition to the cooling water, the fluid delivery system 
also consisted of argon purge gas. Argon was used as a purge gas 
with a flow rate of 25 cm 3 /rain. using a HOKE (Model 133SCAB) 
precision micrometer ing valve. 

For these set of experiments (set *2). a Valley Forge linear 
programmer was used to control power to the preheater. This 
controller ramped the preheater voltage based on a temperature 



furnace arrangement 
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messuremenc fro. . fyp. K thermocouple. The specif Ic.clon for this 
cnecmocouple .r. given In Appendix A The thermocouple was Ins.rtad 
becween Che preheater and Che zlrconla heating element and v.s 
manually removed once the zlrconla Ignition was achieved The 
rlrconla main heater element was connected to a variable 
transformer through a current meter. Th, transformer was calibrated 

using *. volcmettr . 

Th. following s.ction discuss., proc.dur.s us.d to .xamin. 

both f urn.c.s . 

3 . 2 EXPERIMENTAL PROCEDURES 

In this s.ction, th. .xp.rim.nt.l proc.dur.s required to 
verify the performance of th. ARTCOR Model 460-15 Tubular Research 
Furnace are discussed. The procedures were devised to accomplish 
several objectives, Initial experiments w.r. designed to 
characterize the performance of the furnace. Th.s. experiments are 
described in section 3.2.1. The results from thes. experiments gave 
indications of the areas for which improvements could be made. 
Auxiliary equipment had to be integrated into the experimental set 

up to permit further investigations. 

There were three major design improvements that were 
considered in this study. The first was the improvement of the 
furnace components for future us. in microgravity experiments 
discussed in section 3.1.2. The containment of the furnace 
components would be required for microgravity use. It should be 
noted that free convection was us.d to cool th. furnac. casing. 
Since the furnace ends were open to th. atmosphere, on. recognized 
that free convection and radiation losses would play major roles 


the furnace power consumption. Attempts to utilize this furnace for 
materials processing would require correction of these limitations. 

Secondly, the new zirconia heater elements were designed, 
based on oxygen sensor technology. Improvements on these elements 
were concentrated on their ability to carry higher currents loads. 

The wall thickness of the alumina was deduced from the numerical 
studies that were performed. A detailed procedure for the 
construction of these elements is given in section 3.2.3, 

In section 3.2.2 experiments were designed to study the 
temperature and heat transfer characteristics of the research 
furnace. The correlation between the optical pyrometer measurements 
and centerline temperatures were studied and for the two 
experimental setups . 

In the following section a general view of the normal 
operation of the furnace is given. 

3.2.1 Commercial Furnace Operation 

In this section, the normal furnace operation as shown in 
Figures 3.2 is reviewed. 

In this set of experiments, the power supply to the preheater 
and main heater were controlled by a preprogrammed Eurotherm Model 
818P controller. The controller was preprogrammed to perform 
specific functions before the start of the experiment. Initially, 
the controller turned the preheater on at 20% applied power. It 
later ramped the applied power to ^0% at a rate of 1% per minute. 

The preheater was then held at ^0% power for a specified 

time to allow the main heating element to reach its conducting 

temperature. The controller would then switch the power supply from 
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the preheater to the main heating element. The power to the main 
heater was then ramped to a preset value and held for a specified 

time . 

After the given time, the controller would ramp down the 
power to the main heating element at a rate of 1% per minute until 
shut off. Before the next start up of the furnace, the controller 
was initialized and, if desired, reprogrammed to go through the 
routine again. 

As discussed in section 3.1.2 above, the furnace was 
controlled by % power to the zirconia main heating element. The 
temperature measurements for the furnace constitute a crucial role 
in determining its performance in materials processing. Since 
accurate process temperatures are required in various materials 
processing experiments, correlation between the pyrometer readings 
and the centerline temperatures of the furnace needed to be 
studied. Also a correlation between the % power to the element and 
the temperature of the element was an important task in determining 
the performance of this furnace. 

The controller was programmed to power the main heater 
element at 15% to 30% applied power and the temperature responses 
of the element studied. The goal was to partially calibrate the 
furnace temperature as a function of applied power. The 
temperatures measured for different applied power outputs to the 
main heater element are given in Figure 3.1. 

To measure the temperature profile using a pyrometer, a 
special mount to accurately position a high temperature optical 
pyrometer to measure wall temperatures in the zirconia furnace was 
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installed. This mount provided a very precise positioning of the 
pyrometer for consistent temperature measurements. A small viewing 
port was provided by the manufacturer at the side of the furnace 

for direct viewing of the zirconia heating element with the 
pyrome ter . 

A*ial variations in temperature could also be used to 
characterize a furnace. To investigate the axial centerline 
temperature gradient, the furnace was held at applied powers of 20% 
and 30%. The temperature was allowed to stabilize before the 
readings were taken along the centerline of the furnace tube. 

In order to investigate how the furnace responded to gradual 
increases in power input, the furnace was ramped at a rate of 1% 
per minute and the temperature recorded. This experiment was 
performed to evaluate the temperature overshoot which indicates the 
robustness of the controller. 

The commercial furnace system (see Figure 3.2) was equipped 
with an ammeter and voltmeter to measure currents and voltages 
applied to both the preheater and main heater at different power 
settings. This experimental setup was designed to study the actual 
power consumed. The unit was ramped from 50-70% applied power 
following recommended procedures. Additionally, separate 
experiments were conducted at 50, 60, and 70% applied power. 

3.2.2 Research Furnace Operation 

For the second set of experiments, the research furnace was 
equipped with two independent power sources as shown in Figure l.U. 
The preheater was powered by a Valley Forge linear programmer at 
the race of 3°C per minute, to about 800°C. The preheater was then 


held isothermally for about an hour. A variable transformer then 
powered the main heating element through a current meter. When the 
zirconia main heating element ignited, as evidenced by the current 
meter, the power to the preheater was disconnected. 

One may deduce that for set #1, the zirconia element may not 
ignite if the setting for the preheater power and soak time were 
not sufficient to bring the zirconia element to its conduction 
temperature. In set #2, the preheater power was kept on for as long 
as the zirconia element had not ignited (monitored by an ammeter) , 
thus assuring the ignition of the zirconia heating element in each 

experimental run. 

During shutdown, the power to the heating element was reduced 

at a rate of 2% per minute to shut off. 

With the furnace equipped as shown in Figure 3.2, experiments 
were performed to compare the thermocouple (centerline) temperature 
measurements to the optical pyrometer viewing the exterior wall of 
che heating elements. 

3.2.3 Construction of Zirconia Furnace Heating Element* 

During this research, experiments were performed on the two 
zirconia elements available for study. Both elements cracked and 
failed after only limited intermittent use. A possible reason for 
the cracking failure may be the periodic use of the furnace in a 
humid environment. Experimental evidence presented in following 
sections suggested that moisture may have absorbed into the 

elements when the furnace was not in use. 

To study the effects of moisture on the life span of the 
heater elements, a number of elements, in addition to the two 
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obtained from the manufacturer, were required. Due to the high 
costs of elements, experimental zirconia elements, based on oxygen 
sensor technology, were designed and the various components 
ordered. One of these experimental elements was used in controlled 
experiments to determine the initial lifetime. Others will be used 
to determine the effects of moisture. 

The improved elements were designed and manufactured from 
high purity zirconia tubes with 0.00635 m (1/4 inches) platinum 
bands on each end. These zirconia tubes were manufactured by 
McDanel Refractory Company. The zirconia elements were designed 
with inside diameters of 0.0142 m (0.56 inches) and lengths of 
0.1143 m (4.5 inches). This larger inside diameter was selected so 
that different alumina crucibles based on the numerical model 
predictions (discussed in chapter 2) could be studied in an effort 
to design a crucible with an optimum wall thickness. 

Platinum-10% rhodium wire measuring 1000 feet in length with 
a 0.003 inch diameter, was manufactured by the Englehard 
Corporation. This wire was cut into 80 2-ft. strands and twisted to 
form the lead wires for the zirconia heating element. The number of 
strands needed was calculated based on the current loads expected 
during the operation of the furnace. The length of the lead wires 
were established to accommodate the construction of about five 
additional heater elements. The lead wires can be reused to 
construct other elements by cutting off the portion that was 
attached to the spent zirconia heating element, thus resulting in a 

decreased construction cost per element. 

The wires were attached to the zirconia heating element to 


form a current carrying junction between the zirconia heater 
element and the lead wires using a special approach. The junction 
was then cured several times using a Watlow furnace. The element 
was initially cured at 110°C. The temperature was then increased a 
a rate of 5°C per minute to 225°C and held at that temperature for 
another hour. This operation ensured that trapped gases would 
diffuse out of the Junction. The temperature was then raised to 
525°C at the rate of 5°C per minute and held there for about one 
hour. At this stage all the gases and moisture were removed. The 
temperature was raised to 1025°C and held for another hour. The 
power was removed from this well insulated furnace and the system 
allowed to cool undisturbed overnight. The element was installed i 
the furnace cavity and operated following the procedures outlined 
in section 3.2.2. 

In the following chapter the results of the experiments and 
the numerical modeling are discussed. 


CHAPTER 4 


RESULTS AND DISCUSSION 


To verify the performance of the ARTCOR Model 460-15 furnace, 

several experiments had to be performed. The commercial furnace was 

initially analyzed using a number of experiments to determine its 

operating parameters. The furnace components were then analyzed and 

redeveloped as needed to improve the furnace's performance. The 

auxiliary systems were included to adapt the furnace for future use 

in material, processing experiments A major development of this 

study was in the design of the main zirconla heating element which 

was developed based partially on results from the numerical 
modeling. 


Investigations on the heat transfer within the zirconla 
furnace and alumina crucibles were further enhanced by the 
numerical analyses performed in this study. Int.gr.tion of these 
numerical and experimental studies is necessary if the furnace is 
to be used for materials processing experiments. 

The numerical results will be discuss.d in section 4.1. The 
experimental results for both the commercial furnace and the 
research furnace will be discussed in section 4.2. 

4.1 ANALYTICAL AND NUMERICAL RESULTS 
411 011 e Dimensional Analysis 

To evaluate the model, several analyses were performed 
following the works of Vlskanta and Crosh (1962), with fixed 
thermal boundaries, while varying the conduction/radiation 
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parameter N - T 3 +n/v r or . 

o /K • For this analysis, as N approaches 0 one ' 

solves a purely radiation problem, while as N approaches Infinity 

on, solves a purely conductive problem, Figure, 4.1, 4.2 .„ d 4.3 

show the temperature profiles which would be generated within a 

thlch alumina crucible at various value of N and fired temperature 

differences. For this analysis only the crucible wall was 

considered. Figure 4.1 shows the results for an extreme tempera- 

sure difference of 1713K (1440«C>. This figure shows that at large 

temperature differences, materials with very high semitransparent 

properties would be required to distort the thermal profile within 

the crucible. Basically, this analysis Is representative of a 

radiation problem with constant emisslvlty followed by conduction 

to the ocher boundary. At lower cemperature differences. Figures 

‘■2 and 4.3, radiation effects within the material begin to take 

affect and only small changes in the materials conduction/radiation 

parameter produce large gradients within the sample. 

The above analysis was undertaken to verify literature 
results, as It would require actual variations In the crucible 
material of construction to vary the conductlon-radlation 
parameter. The alumina crucible modeled her. Is a component of an 
actual furnace and Is mad. of high purity alumina. Ac the furnace 
operating temperature, the crucible has almost a constant value for 
N The numerical problem addressed In this study was principally 
concerned with cemperature measurements of an alumina crucible 
using an optical pyrometer and prediction of the thermal profile 
Within a liquid iron sample based on these measurements. Refer to 
f-gJte 2.2 to review the system being analyzed. Optical pyromecers 
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generally require an average emissivity value to calculate surface 
temperatures, and thus an opaque boundary is assumed. In Figure 
4.4, is shown a comparison of an opaque alumina crucible with an 
average emissivity with heat transferred within the system purely 
by conduction, to the model developed in this study for a 
se mi transparent material with the heat flux matched at the 
crucible/liquid metal interface. In Figure 4.4, this interface is 
located at radial distance of 1 mm. Since an optical pyrometer 
assumes an opaque boundary, it would over estimate the wall 
temperature by about 5°C. If one used the opaque model for open 
loop control, they would calculate a sample temperature 2-3°C lower 
than that of the semitransparent model. The semitransparent model 
(Figure 4.4) shows the expected hot spot within the crucible wall. 
This hot spot results from the various depths to which given 
frequencies of radiation can penetrate before being absorbed by the 
material. The degree to which the profile within the crucible 
influences the gradient within the liquid iron is a function of 
crucible thickness. Here a 1 mm crucible was assumed. Figure 4.4 
suggest a possible solution for high temperature furnace control 
problems, ones in which direct temperature measurement techniques 
are impossible to employ due to material constraints. A combined 
approach would correct the optical pyrometer's temperature 
measurement using a semitransparent model to predict, in real time, 
the temperature profile within the sample. Thus, the power to the 
furnace element would be adjusted after a computer calculation of 
the predicted temperature versus the set point temperature was 
made. It should be noted that opaque boundary is commonly used to 
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Figure 4,4 Temperature profiles for the Alumina Crucible and 
Liquid Iron Sample. — represents the profile 
for the semitransparent case while — — •represents 
the profile for the opaque assumption. 
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solve combined conduction/radiation problem# in computer programs 
such as SINDA (Smith 1971) . The one dimensional model failed to 
predict the thermal profiles within crucible walls as a function of 
crucible thicknesses, therefore a two dimensional model w*s 
required . 

Vi. 2 Tvo Dimensional Analysis 

One of the main objectives of the two-dimensional analysis 
was to examine the effect of the variation in temperature profiles 
within crucible walls of varying thickness, with the heating 
element at a constant temperature of 1600°C. The results are shown 
in Figure 4.5 for different values on N (the radiation-conduction 
parameter) . Here N is proportional to the ratio of crucible to 
liquid iron thicknesses, namely 1:3, 1:1, and 4:1. Therefore the 
variations in temperature profiles are a result of the total heat 
capacities of the crucibles and the heat dissipation in the axial 
direction within the liquid iron sample. Since heat is dissipated 
from both the top and bottom of the sample, the temperature 
profiles towards the ends will be more pronounced chan the profiles 
within the center of the model. 

When the relative thicknesses for the alumina crucible and 
the liquid iron samples are the same (see Figure 4.6) it was noti- 
ced that conduction in the liquid iron is strongly influenced by 
che thermal profile developed within the crucible wall. Notice chat 
as che crucible wall gets thinner, the thermal profile increases to 
a maximum almost at che liquid iron interface. This suggests chat 
an optical pyrometer observing che crucible wall would be largely 
measuring the temperature of che liquid iron surface. 
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In Figure A. 7, che ratio of the crucible wall thickness to 
that of the iron sample is 1:3. The temperature at the 
crucible/iron interface is much greater than for the previous 
cases. This is desirable since for a given power to the furnace, 
one would like to achieve a high temperature at the liquid iron 
interface while ignoring the crucible wall. This means that the 
radiative effects within the crucible are minimized with respect to 
che bulk sample. 

4.2 EXPERIMENTAL OBSERVATIONS 

In section A. 2.1 the results from che experiments performed 
and the problems encountered while using the commercial zirconia 
furnace are discussed. The experiments performed using the research 
furnace, including the design improvements are discussed In section 
A. 2. 2. 

4.2.1 Commercial Furnace 

The initial experiments were designed to study the heat 
transfer characteristics of the commercial furnace. Figure A. 8 
shows the centerline temperatures measured 1 Inch from the bottom 
of the furnace by a thermocouple. These temperatures were matched 
to the pyrometer readings by increasing the Input power to the 
furnace. The results are given in Figure A. 9. The results show a 
factor of two for the power requirements. 

Figure A. 10 shows the measured axial temperature profile for 
20% applied power, while Figure A. 11 shows the results at 30% 
applied power. Axial gradients of approximately 523K (250°C) per 
inch were determined for both figures. Temperatures were found to 
be higher at the furnace cop than at equivalent distances from the 
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Figure 4.9 Temperature as a function of Applied Power for the 
ARTCOR Model 460-15 Furnace as Measured by a Uater 
Cooled OS-IOOO-HT Optical Pyrometer. 
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Figure U . 10 Axial Temperature Profile for the ARTCOR Model 
^*60-15 Furnace at 20% Applied Power on the 
Controller . 





furnace bottom as would be expected for a freely convecting medium. 

Figure 4.12 shows the temperature response to a gradual 
increase in power at the rate of 1% per minute. It was found that 
the ratio of the increase in temperature to the applied power was 
approximately 1:15. A temperature overshoot of about 473 AK (200°AC) 
was recorded. 

After numerous experimental runs of the ARTCOR Model 460-15 
Research Furnace, several problems were encountered. When the 
furnace was left for 48 hours or more, the main element would not 
Ignite using the normal procedure of preheater setting of 40% 
applied power for 0.6 hours. It may be that the main heater element 
had a tendency to adsorb water vapor from the atmosphere, thus 
requiring a ^higher input of power by the preheater in successive 
runs . 

The commercial furnace (shown in Figure 3.2) would not 
operate above a 50% power settings using the two commercial 
elements. The current and voltage achieved at 50% applied power 
were 3.7 amps and 80.0 volts, respectively. Identical current and 
voltage results were obtained at 60 and 70% applied power on the 
controller. Thus, a maximum power had been achieved for the 
commercial element . 

Two commercial elements were used during a portion of this 
study. Both the elements failed prematurely compared to the 
manufacturer's suggested lifetime. The first zirconia element 
lasted for a total of 15 hours at main heater power settings 
varying from 15% to 50% applied power. The power setting of 50% 
corresponds to about 4 amps and 80 volts across the element, or 320 
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waccs. The second element was initially started using a 43% 
preheater power for 0.8 hour. This procedure was able to start the 
aain heater element, but may have been the cause of its shorter 
lifetime as will be explained later. This element lasted for a 
total of on.y 3 hours at main element powers of 40% to 50%. 

Figure 4.13 is a photograph of the two zirconia main heater 
elements. Element 1 is the shorter of the two elements shown. Under 
a microscope, both elements showed defects that appeared as glassy 
regions with a channeled structure at the points where the cracks 
could have initiated. Both cracks seemed to have initiated in the 
hottest zone of the furnace. 


A possible explanation for the difference in lifetime between 
the two elements studied, may be that the first element received a 
slow cure while the second received a much faster cure. Slow curing 
may have removed water from the ceramic and further sintered the 
ceramic, providing greater strength. Additionally, preheating at 
1073 - 1173K (800 - 900°C) may have partially annealed the ceramic 
thereby reducing stresses within the element. 


Another possible explanation could be that the element, 
relative to its final operating condition, may be partially green. 
The elements may contain trace water within the matrix. Element #2 
received both a rapid cure and a higher applied power than element 
*1. Perhaps element #1 had sufficient time for the water to diffuse 
through the nonporous zirconia matrix prior to heating at a lower 
power setting. 


Pure water may act as an insulator if trapped within the 
zirconia element. Water trapped within the matrix at the high 
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igure 4.13 Zirconia Main Heating Elements. Element #1 
snorter of the two elements which lasted a 
of 15 hours, while element //2 lasted only 


is the 
total 
3 hours. 



operating temperatures of the zirconia main heating element would 
convert to steam within the matrix resulting In an extremely high 
local pressure gradient. Element #2 may have been manufactured at 
the same time as element «1 and was stored locally for about 8 
months. During this time it may have accumulated additional 
moisture from the atmosphere. Water is again suggested in the extra 
power required to heat the elements after sitting in a humid 
environment. This is partially supported by the microscopic 
investigation discussed above. 

If we assume that water adsorbed into the element, it needed 
to be evaporated prior to ignition. This may have caused a pressure 
gradient within the element which resulted in a high stress 
distribution that lead to cracks and ultimate failure of the 
element. 

After the preliminary studies were performed using the 
commercial heating elements and the results analyzed as above, the 

study focused on efforts to design improved zirconia elements and 
furnace component parts . 

4.2.2 Research Zirconia Furnace 

In order to analyze why the element's lifetime was shorter 
than expected, several controlled experiments needed to be 
performed; therefore a more economical design of the furnace 
elements for research was required. 

In an effort to reduce the cost of replacement elements and 
to study the elements in greater detail, preliminary studies were 
performed on techniques to design new experimental elements, that 
reduced the cost and performed as required. There were two basic 
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co„po„ e „ ts bp th. .am heating element. These were ba.lc.lly 2 high 
temperature high current lead vires and a high density high purity 
piece zirconia tubing. 

The new main heater elements, were manufactured from high 
Parity rlrconla tubes with 1/a inch platinum bands on each end. 
Platinum bands were added to permit better contact between th, 
tlrconla tube and the platinum- rhodium lead wires. The method of 
attachment of th, wires to th, platinum band was devised based on 
oxygen sensors technology. The lead wires, composed of 90% platinum 
- 10% rhodium alloy, were selected to carry th. required current at 
the maximum operating conditions Th, Pt band/past, combination 
provided . good electrical connection to th. rlrconla element. 

The experimental runs were designed to study several aspects 
-he new zirconia heating elements. These newly developed 
elements (with layout shown In Figure 3.4) were powered to values 
of 70% power at currents of 6.0 amps and voltages of 80 volts. In 
-r.-s set or ex-^rjMr.ts (set *2), the furos.ce equipped with tw< 

independent power supplies, one to power the preheater and the 
other to power the main heater element. This configuration 
eliminated the uncertainties in igniting the main zirconia heater 
element, that were caused by the element not reaching its 
conduction temperature (1073K) . 

Figures 4.14 and 4.15 show the results of one of these 
experiment. In Figure 4.14 is shown the current history from the 
element operated at a constant voltage of 98 VAC. Notice that the 
current rises as a function of time due to the fact that as the 
element heats its resistance continually decreases (Figure 4.15) 
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igure 4.14 Variation of Current with Time for a Constant 
Voltage of 98 VAC. 
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Figure 



4.15 Variacion of Resistance with Time for a Constant 
Voltage of 98 VAC. 
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from about 445 ohms to a final resistance of 21.2 ohms. This 
corresponds to a current of 4.62 amps after about 1 1/2 hours of 
operation. The new zirconia element operated for over 25 hours 
without failure. The commercial elements carried a measured maximum 
current of 4.6 amps, while the new element was found to carry 
currents of over 7 amperes. Higher current loads measured in the 
new element may resulted from techniques used to attach the lead 
wires . 

The centerline temperatures that were measured for the the 
Improved zirconia heating elements are shown in Figure 4.16. Also 
shown in Figure 4.16 is a comparison of the temperature readings 
between the thermocouple measurements of the centerline 
temperatures within the alumina crucibles and the pyrometer 
readings through the viewing hole. The reason for this wide 
difference in the readings was due to differences between the 
aperture of the pyrometer and the diameter of the viewing hole 
provided by the manufacturer. This difference would be very 
significant if the furnace in its current configuration were to be 
used in materials processing experiments. 


TEMPERATURE 




CHAPTER 5 


CONCLUSIONS AND RECOMMENDATIONS 


5.1 CONCLUSIONS 

In the one dimensional radiation/conduction analyses it was 
shown that the transmissivity of a material has a significant 
effect on the temperature profile when radiation is considered at 
such high temperatures. Considering the physical model, it was 
shown that the crucible wall appears thermally thinner than its 
physical thickness. Since the hot zone is within the crucible wall, 
temperature measurements made using optical pyroraetry need to be 
adjusted if adequate control is to be achieved. It has been 
demonstrated that the variable radiation properties are significant 
in the analysis of radiation/conduction problems at such extreme 
temperatures . 

The two dimensional model proved to be very successful in 
modeling the furnace crucible. The results have provided a greater 
detailed understanding of the components of this advanced high 
temperature furnace. The parametric study of the ratio of the 
crucible wall to the liquid Iron sample, has shown that there is a 
significant Increase In the temperature at the interface (between 
the liquid iron sample and the crucible wall) with decreasing 
crucible thickness. Also, the radiation effects on the sample are 
greatly reduced for thin crucible walls. This allows a greater 
portion of the iron sample to attain a more uniform radial 
temperature profile relative to the centerline. For a thick wall, 
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an optical pyrometer would have measured the average temperature of 
the crucible wall. Optical pyrometer measurements for a thin 
crucible wall would approximate the interfacial temperature at the 
liquid iron sample. The general conclusion from this portion of the 
study is to select as thin a crucible as will tolerate the loads 
applied during processing. 

Initial performance data for a new zirconia element design 
was presented, and results show consistent performance and greater 
power handling capabilities. As discussed in the results, the 
higher current loads are a result of the improved junction used for 
this zirconia element. Procedural modification followed for the new 
element, indicate that the preheater power should not be shut off 
until the zirconia element has ignited. Therefore one has to either 
have two independent power supply sources for the preheater and 
main heater, or have a feedback control scheme that will monitor 
the current in the main zirconia element prior to disconnection of 
the preheater power. 

Centerline temperature measurements using thermocouples 
recorded temperatures twice as high as the pyrometer measurements 
for equal power inputs. After modification of either the pyrometer 
focal length or increasing the diameter of the viewing hole, better 
correlations between these two temperature measurements will be 
possible. With the integration of the numerical modeling and the 
temperature measurements from the optical pyrometer, it may be 
possible to develop a feedback control scheme for use with these 
furnaces. Electronics with variable eraissivities accessible from a 
digital computer would greatly aid in the accurate measurement of 
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Che wall temperature when using commercially available optical 
pyrometers . 

5.2 RECOMMENDATIONS 

Several improvements could be made in the area of control for 
this type furnace. Detailed studies are needed to verify the extent 
and sensitivity of the zirconia elements to water vapor. 

Experiments are needed that will measure the centerline temperature 
of the furnace while solidifying several pure metal systems. This 
approach will provide a consistent reference for numerical modeling 
and to both calibrate the experimental apparatus. 

A suggested experimental setup is shown in Figure 5.1. For 
this open loop control scheme, the computer would contain a 
numerical model capable of predicting the temperature profile of 
the furnace sample based on accurate optical pyrometer 
measurements . The computer could also estimate the average 
emissivity based on the spectral bandpass of the emitting wall as a 
function of temperature. This value would serve as input to the 
optic pyrometer, thus minimizing another source of potential error 
that occurs when an average emissivity is used over the entire 
range of operating temperatures. In addition the computer should be 
equipped with current monitors to control both the preheater and 
main heater power supplies. 



Figure 5.1 Recommended Layout of the Experiment including 
Auxiliary Systems . 
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APPENDIX A 


TEMPERATURE MEASUREMENT DEVICES 

I. THERMOCOUPLES 
TYPE S 

Platinum- 10% Rhodium • Platinum 
Temperature Range 335 * 1750°C 
Error 1 2 . 3°C 
Resolution 1°C 

TYPE K 

Chromel * Aluael 
Temperature Range 200 - 1250°C 
Error 1 2.0°C 
Resolution 1°C 

II. OPTICAL PYROMETER 

TYPE OS - 1000 -HT 
Optics 

Spectral Response 2.0 • 2.5 microns 
Resolution 50 : 1 ( 1°) , 90% energy 

Minimum Spot Size 1.5 cm at 76 cm 
Accuracy 

-1% of reading, -1 digit at 25°C head temperature 
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APPENDIX B : FORTRAN PROGRAM FOR ONE DIMENSIONAL 
RADIATION/CONDUCTION MODEL 

PARAMETER (ND - 100, ND2 - 200) 

EXTERNAL AFUNC, THERM! , THERM2 , ADS IMP 
DIMENSION A(ND.ND) ,B(ND) ,T(ND) ,V(ND) , 

$ DIFF(ND) , AFT(8) , ER ( 2 ) , TC (ND) , QT (ND) ’ 

$ WW(ND) , EE(ND) , FT(ND2 ) , DIST(ND) 

C 

SIGMA - 5 . 669E-12 
TOL - 1.0E-6 
IER - 0 

OPEN (6, FILE - 'CRUCL.DAT*, STATUS - 'NEW') 

C INITIALIZE THE MATRICES 
DO 1 I - 1, ND 
B(I) - 0.0 
T( I) - 0.0 
W(I) - 0.0 

wv(i) - o.o 

DO 1 J - 1, ND 
A(I,J) - 0.0 
1 CONTINUE 

C RANGE OF WAVE LENGTHS AND INITIAL GUESS TEMPERATURES 
WRITE(* , *) 'ENTER T(l) T(ND) ' 

READ(* , *) T ( 1 ) , T(ND) 

WRITE(* , *) 'ENTER SL AND BL ' 

READ(* , *) SL, BL 
H - 1/ND 
DIST(l) - 0.0 
DO 6 I - 2, ND 

DIST(I) - DIST(I-l) + H 
6 CONTINUE 

COUNT - 0.0 

200 IF (COUNT. GT. 0.0) THEN 

IF (FLX2.GT.FLX1) THEN 
T(ND) - T(ND) + 3.0 
P - T(ND) 

ELSE 

T(ND) - T(ND) * 2.5 
P - T(ND) 

END IF 
ELSE 

WRITE ( 6,*) 'COUNTER - '.COUNT 
END IF 

DO 5 I - 1, 8 

AFT(I) - ADS IMP (AFUNC , SL , BL, TOL, IER , I) 

5 CONTINUE 

CALL RAD(ND , T , B , H , SL , BL, ND , FLX1) 

WRITE(*,*) 'FROM RADIATION FLUX AT JOINT - ' FLX1 
DO 2 I - 1, ND ' 

FT(I) - T(I) 

2 CONTINUE 

CALL CONDUCT ( ND , T , H , HH , ND , FLX1 , FLX2 ) 
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WRITE(*,*) 'FROM CONDUCTION FLUX AT JOINT - ' FLX2 
DO 3 I - 1, ND2 
FT(I+10) - T( I ) 

3 CONTINUE 

COUNT - COUNT + 1 
IF(COUNT . GT . 100 , 0)G0TO 300 
C CONDITION FOR STOP 

I F( COUNT . GT . 0 . 0) THEN 
ER(2) - ER( I) 

ER(1) - ABS ( FLX2 - FLX1) 

ERA - ABS(ER(1) • ER(2)) 

IF (ERA.CT. 1.0E-2) THEN 

WRITE (*,*)' ERROR IN FLUXES '.ERA 
GOTO 200 
ELSE 

URITE(*,*) 'ERROR AND FLUX1 FINAL '.ERA FLX1 
END IF 
ELSE 
END IF 

300 WRITE(* , *) '# OF ITERATIONS - '.COUNT 

WRITE(*,*) 'ERROR IN FLUX AT JOINT - ' , ERA 
WRITE(*,*) ' FLUX AT RAD/COND JOINT - ' ' FLX1 FLX2 
DO 4 I - I, ND2 

WRITE(* , *) I , ' - ' , FT(I) 

4 CONTINUE 
STOP 
END 

C 

SUBROUTINE RAD ( ND . T , B , H , SL , BL , ND , FLX1 ) 

DIMENSION A(ND , ND) , B(ND) , T(ND) , W(ND) , DIFF(ND) , AFT(8) 
$ TC(ND) ,QT(ND) ,VU(ND) , EE(ND) , FLXR(ND) 

$ PT(ND) 

T(l) - Q 
T(ND) - P 
PT(1) - T( 1 ) 

COUNT2 - 0.0 
SIGMA - 5.669E-12 
TOL - 1.0E-6 
IER - 0 
WV(1) - 0.0 
WWW - 0.0 
DO 423 I - 2 , ND 

WW(I) - H + WWW 
WWW - WWW +H 
423 CONTINUE 
C 

Z - (T(l) - T(ND) )/ND 
DO 4 I - 2, ND-1 
T(I) - T( 1) • I*Z 
PT(I) - T(I) 

4 CONTINUE 

DO 41 I - l.ND 
QT( I ) - T(I)/T(1) 

41 CONTINUE 


c 


MM - 1 

adjusting guess values 
IF(MM . GT . 0)GOTO 212 
37 DO ^2 I - 2, ND-1 

I F ( PT ( I ) .GT.T(I) ) THEN 


PT(I) - T ( I ) 
T( I ) - T( I) + 


0. 5*DIFF(I) 


QT( I ) - T(I)/T(1) 
ELSEIF(PT(I) . LT . T ( I ) ) THEN 


PT ( I ) - T ( I ) 

T ( I ) - T(I) - 0 . 5*DIFF ( I ) 
QT ( I ) - T(I)/T(1) 


U2 

212 


22 


ELSE 
END IF 
CONTINUE 
Y-1.0 

DO 22 I - l.ND 

TC ( I ) - 0 . 01*THERMl (T ( I ) ) 

CONTINUE 
C GENERATE MATRIX [A] and [B] 

A(l.l) - l-° 

A(ND.ND-l) - TC(ND) + TC(ND-l) 

A(ND,ND) - -A(ND.ND-l) 

D °A(l!l’l)'- 1 '<O l 5*(TC(I-l) + TC( I) ) )/(H**2) 
A( I , I ) -- (0 . 5*TC(I- 1) + TC(I) + 

$ 0 . 5*TC(I+1) )/ (H**2) 

A( I , 1+1) - (0 . 5*(TC(I) + TC(I+1) ) )/(H**2) 

2 CONTINUE 


CALL FORCE ( ND , QT , W , B , AFT , S L , BL , TOL , I ER ) 
B(l) - 1.0 

B(ND) - AFT(2)*SIGMA*(T(ND)**M 
CNTERM -( SIGMA* (T(l)**3) )/AFT(3) 

DO 5 I - 2, ND-1 

B(I) - -1 .0*(CNTERM*B(I) ) 

5 CONTINUE 


C 

CALL GAUSS (ND , A , B) 


C 


666 


8 


13 


DO 666 I - 1, ND 
QT(I) - B(I) 

T ( I ) - T(1)*B(I) 

CONTINUE 
DO 8 I - 1, ND 

W(I) - T(1)*B(I) 

WRITE (* , *) WW( I ) , ' . ' ,W(I) 


T(ND- 2) ) / ( 2 . *H) 


CONTINUE 

$ THERM1(T(ND)) 

DO 13 I- 2 , ND 

DIFF(I) - ABS ( PT ( I ) - T(D) 
CONTINUE 



C CONDITION FOR STOP 
MM - MM+1 

COUNT 2 - COUNT 2 +1.0 
IF (COUNT2 . GT . 20 . 0) GOTO 300 
DO 14 I - 1, ND 

IF(DIFF(I) .GT.0.5) COTO 37 
14 CONTINUE 

300 WRITE(*,*)' TOTAL ITRS FOR RADIATION ' . COUNT2 

RETURN 

END 

C 

SUBROUTINE CONDUCT (ND , T , H , HH , ND , FLX1 , FLX2 ) 
DIMENSION A(ND , ND) , B(ND) , T(ND) , U(ND) , DIFF(ND) , 

$ BFLX(2) ,TC(ND) ,QT(ND) ,UU(ND) , EE(ND) ’ 

$ FLXC(ND) 

C 

T(l) - T(ND) 

COUNT 3 - 0.0 
SIGMA - 5 . 669E-12 
TOL - 1.0E-6 
IER - 0 
UU(1) - 0.0 
UUU - 0.0 
DO 423 I - 2 , ND 
UU(I) - HH + UUU 
UUU - UUU +HH 
423 CONTINUE 
Z - 1.0 

DO 4 I - 2, ND 

T(I) - T(I-l) - Z 
4 CONTINUE 

C GENERATE MATRIX [A] and [B] 

DO 41 I - l.ND 
QT(I) - T(I)/T(1) 

41 CONTINUE 
MM - 1 

C ADJUSTING GUESS VALUES 
NDD - ND - 1 
I F ( MM . GT . 0 ) GOTO 212 
37 DO 42 I - 2, ND- 1 

IF(BFLX( 1) . GT . BFLX(2) )THEN 

IF( I . EQ . NDD)T(ND) - T(ND) ' - 0.5 
T ( I ) - T( I ) - 0 . 1*DIFF ( I ) *T ( 1 ) 

EE(I) - QT( I ) 

QT(I) - T(I)/T(1) 

ELSE 

I F ( I . EQ.NDD)T(ND) - T(ND) + 0.5 
T ( I ) - T( I ) + 0 . 1+DIFF ( I ) *T ( 1 ) 

EE(I) - QT( I ) 

QT( I ) - T(I)/T(1) 

ENDIF 

42 CONTINUE 

212 Y-1.0 
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DO 22 I - l.ND 

TC ( I ) - 0 . 01*THERM2 (T(I) ) 

22 CONTINUE 

A(l,l) - 1.0 

A(ND.ND-l) -- (TC(ND)+TC(ND* 1) ) 

A(ND.ND) - *A(ND , ND* 1 ) 

A(ND.ND-l) - 1.0 
DO 2 I - 2. ND-1 

A(I.I-l) - (0 . 5*(TC( I - 1) + TC(I) ) 

A(I.I) - - (0. 5*TC( I- 1) + TC ( I ) + 0 . 5*TC( 1 + 1) ) 

A( I , 1 + 1 ) - (0 . 5*(TC ( I ) + TC ( 1+1 ) ) 

2 CONTINUE 

B(l) - 1.0 
DO 5 I - 2 , ND 
B(I) - 0.0 
5 CONTINUE 

C 

CALL GAUSS (ND , A , B) 

C 

DO 666 I - 1, ND 
QT( I ) - 8(1) 

T(I) - B(I)*T(1) 

666 CONTINUE 

DO 8 I - 1, ND 

U(I) - T ( 1 ) *QT ( I ) 

WRITE (+ , *)UU( I ) , ' , ' ,U(I) 

8 CONTINUE 

BFLX(2)-BFLX(1) 

FLXC(l)- ( • 3 . *T(l) + 4 . *T ( 2 ) - T ( 3 ) ) / ( 2 . *HH ) 
WRITE(*,*) 'CONDUCTION flux at 1 - -,FLXC(1) 
FLXC(ND)- ( 3 . *T(ND) • 4.*T(ND-1) + 

$ T(ND-2))/(2.*HH) 

WRITE (*,*) 'CONDUCTION flux at ND - ' , FLXC(ND) 
FLX2 - FLXC(l) 

BFLX(l) - FLXC(ND) 

FLX2 - THERM2(T(1) )+FLXC(l) 

DO 13 I- 2 , ND 

DIFF(I) - ABS(QT(I) - EE(I)) 

13 CONTINUE 

C CONDITION FOR STOP 

WRITE(*,*) 'CONDUCTION Iterations - ' ,MM 
MM - MM+1 

COUNT 3 - COUNT 3 +1.0 
IF (COUNT3 . GT . 20 . 0) GOTO 300 
DO 14 I - 1, ND 

IF(DIFF(ND) .GT.1.0E-6) GOTO 37 
14 CONTINUE 

300 WRITE (*,*) 'TOTAL ITRS FOR CONDUCTION — ' COUNT 3 
RETURN 
END 
C 

SUBROUTINE FORCE (ND , ST , W , SS , AFT , SL , BL , TOL , IER) 
DIMENSION ST(ND) ,W(ND) t SS(ND) ,AFT(8) 

REAL Y , XY , AFUNC 


SIGMA - 5.669E-I2 
DO I I - 1, ND 
C 

SSE - 0.0 
C 

C GAUSSIAN QUADRATURE 
C 


C 

C 


C 


2 

I 


C 


IF (J.EQ.l) THEN 

Y - -0.57735027 
LAMDA - SL 

ELSE 

Y - 0.57735027 

LAMDA - BL 
END IF 


XY - 5 . 0E-5*Y + 5.0E-4 

A - AFT (1)*(ST(I) **4 ) 

B - ((1.0 - AFT ( 2 ) ) *0 . 1281 )/ 

$ (2 . 0*0 .128 1* (1.0 • AFT ( 2 ) ) 

$ + ( (AFT(l) )*SIGMA*(ST(I)**4) 

$ *( (ST( 1 ) *1600) **4) ) ) 

C - AFT(3)*(ST(I)**4) 

D - ( AFUNC (1.0,8))*(ST(I) **4 ) * 

$ ((AFUNCQCY,1))-AFT(4)) 

F - 2 . 0* ( 1 . 0 - (AFUNC (XY , 6 ) ) ) 

G - (AFUNC (XY, 6) )*(ST(I)**4) 

H - 1.0/0.1281 


AILOG 
SSE - 

$ 

$ 

$ 

$ 

$ 

CONTINUE 
SS(I) . 
CONTINUE 
RETURN 
END 


- 1.0 - XY 

SSE + 5.0E-5 * (4.0 * G 

- 2.0 * (AFUNC (XY , 2 ) ) * 

(G + 2 . 0*F*B* ( (-G - F*A) - H* (C-A) ) + A) 
+ 4 . 0*B* ( -G - F*A - H*(C-A) ) 

* ( 1.0 + (0.5772 - 1.0 + alog(AILOG) ) * 
(l.-XY) - ((1. -XY)**2)/2.0) - 2*D ) 

- SSE 


REAL FUNCTION AFUNC(X.ID) 

REAL X 
INTEGER ID 

5 GOTO (10,20,30,40,50,60,70,80,90) ID 

10 AFUNC - -0.5772 - alog(X) + X - (X**2)/4.0 
goto 90 

20 AFUNC - 1.0 + (0.5772 - 1.0 + alog(X))*X - 
$ (X**2)/2.0 

goto 90 

30 AFUNC - 0.5 • X + 0.5*(-0.5772 + 1.5 - alog(X) )*X**2 
$ + (X**3)/6 .0 
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goto 90 

AFUNC - ((-0.0001 * (X**5) + 0 . 0031*(X**4) - 
$ 0 . 052 3* (X**3 ) + 0.3834 * ( X**2 ) - .9968 * X + 

$ .8177) * (1.003**4) * (-0.5772 - alog( 1 . 003E- 5) 

$ + 1.003E-5 - ( 1 . 003E- 5**2) /4 . 0) 

$ - ( -0 . 0001*(X**5 ) + 0.0031 * (X**4) - 0.0523 * 

$ (X**3 ) + 0.3834 * (X**2 ) - .9968 * x + .8177) * 

$ (1.003**4) *(-0.5772 - alog(X) + X - (X**2)/4.0)) 
goto 90 

50 AFUNC - (1.0 + (0.5772 - 1.0 + alog(X))*X - 

$ (X**2)/2.0) * ( - 0 . 0001* (X**5 ) + 0 . 0031*(X**4) - 

$ 0.0523*(X**3)+ 0 . 3 8 34* ( X** 2 ) - ,9968*X + .8177) 

goto 90 

60 AFUNC - - 0 . 0001* (X**5 ) + 0 . 0031*(X**4) - 
$ 0 . 0523*(X**3) +* 

$ 0. 3834*(X**2) • , 9968*X + .8177 

goto 90 

70 AFUNC - . 0004*X + .4432 
goto 90 

80 AFUNC - 1. 0/(0, 5 - 1.0 + 0. 5*(-0. 5772 + 1.5 - 
$ alog( 1 . 0) )*1 . 0**2 + ( 1 . 0**3)/6 . 0) 

C 

90 RETURN 
END 
C 

C THERMAL CONDUCTIVITY OF ALUMINA 
C 

REAL FUNCTION THERMl(X) 

REAL X 

THERM 1 - . 0019*X +6.0 
RETURN 
END 
C 

C THERMAL CONDUCTIVITY OF IRON SAMPLE 

C 

REAL FUNCTION THERM2 (X) 

REAL X 

THERM2 - . 0219*X +2.6 

RETURN 

END 

C 

C THIS FUNCTION CALCULATES THE INTEGRAL OF THE FUNCTION IN THE 
C REAL FUNCTION 'AFUNC', USING ADAPTIVE QUADRATURE WITH SIMPSON'S 
C RULE . 

C A - LOWER BOUND OF THE INTEGRAL 

C B - UPPPER BOUND OF THE INTEGRAL 

C TOL - MAXIMUM ERROR OF THE QUADRATURE 

C IER - ERROR. IER - 129 IF THE ALGORITHM COULD NOT CONVERGE 

C 

c 

REAL FUNCTION ADSIMP(AFUNC , AA, BB, TOLL, IER.Q) 

C 

REAL AA, BB, TOLL, AFUNC 
INTEGER Q 



DOUBLE PRECISION EST , R( 20 ) , V ( 20) , FR( 20) , 

* FS(20) , FT(20) , FU(20) , FV( 20) , ERR , S2F, SF, HMIN TOL 

LOGICAL FIN ' ’ 1UL 

INTEGER LEV.IER 
C INITIALIZE 


IER - 0 
A - DBLE(AA) 

B - DBLE(BB) 

TOL - DBLE ( TOLL) 

HMIN - (BB- AA)/2 . D0**19 
EST - 0.D0 


R(l) - A 
V(l) - B 
IF (A.EQ.B) 
LEV - L 
FR(1) 

FS(1) 

FT( 1) 

FU(1) 

FV(I) 

FIN - 


RETURN 


DBLE(AFUNC( SNGL( 
DBLE(AFUNC( SNGL( 
DBLE(AFUNC( SNGL( 
DBLE ( AFUNC ( SNGL( 


DBLE(AFUNC(SNGL(V(1) ) ,Q) ) 
.FALSE. 


R(l) ) ,Q) ) 

(V( 1) -R( 1 ) )/4 . DO + R ( 1 ) 

( V ( 1 ) -R(l) )/2 .DO + R ( I ) 
(V(l) -R(l) )*. 75D0 + R( 1) 


C 

C ENTER THE MAIN LOOP 
C 

100 IF (FIN) RETURN 
C 

C CALCULATE S2(F) AND SF(F) 


■ Q) ) 

■ Q) ) 
.Q) ) 


C 

C 

C 

C 

C 

C 


S2F - (V(LEV) -R(LEV) ) / 1 2 . D0*(FR(LEV)+FV(LEV)+2 . D0*FT(LEV) + 
* 4.D0*(FS(LEV) + FU(LEV) ) ) 

SF - (V(LEV) -R(LEV) )/6 . DO*(FR(LEV)+FV(LEV) + 4 . D0*FT(LEV) ) 

GET THE TOLERANCE FOR THIS STAGE 

ERR - 2 . DO*TOL*(V(LEV) -R(LEV) )/(B-A) 

IF THE ERROR IS TOO HIGH, GO TO THE NEXT SMALLER STAGE 
OTHERWISE, DO THIS STAGE. 


IF(DABS(S2F-SF) .GT. ERR) GO TO 200 
EST - EST ♦ S2F 
ADSIMP - SNCL(EST) 

C 

C CHECK FOR ENDING CRITERIA. IF LEVEL - I AND V(LEV) - B, STOP 

c 

IF ( ( LEV . GT . 1 ) . OR . ( V ( LEV ) . LT . B ) ) GO TO 150 
FIN - .TRUE. 

GO TO 100 
150 CONTINUE 
C 

C OTHERWISE. SET UP FOR THE NEXT STACE 
C 

R(LEV - 1) - V(LEV) 

FR(LEV - 1) - FV(LEV) 

FT (LEV - 1) - FU(LEV - 1) 



FS C LEV • 1) - DBLE ( AFUNC ( SNCL( ( V( LEV - 1 ) • R( LEV - 1 ) ■> 

>u D0+R(LE7- 1 ^ O'- 

FU ( LEV • i ) - DBLE ( AFUNC < SNGL( ( V ( LEV- 1 ) • R( LEV • 1 ■ ; * *:00 
- R ( LEV • 1 .. 

LEV - LEV - 1 
GO TO 100 

: IF IT DIDN'T WORK OUT . * GO TO THE NEXT HIGHER LEVEL 


0 SET UP THE BOUNDARIES AND FUNCTIONS FOR THE NEX1 


VEL CALCULATE 


IF (LEV. EQ. 20) GOTO 225 
R( LEV+1 ) - R ( LEV ) 

V( LEV+1) - ( R ( LEV ) +V ( LEV ) ) / 2 . DO 
FR ( LEV+1 ) - FR ( LEV ) 

FT ( LEV+1 ) - FS ( LEV ) 

FV(LEV+1 ) - FT (LEV) 

FS ( LEV+1 ) - DBLE ( AFUNC ( SNGL( (V( LEV+1 )- R( LEV+1 )) /u DO - 
* R( LEV+1) . Q) 

FU ( LEV+1 ) - DBLE (AFUNC ( SNCL( ( V ( LEV+1 ) - R( LEV+I ) ) * . 7 5D0 - 
- R ( LEV + 1 ) , Q ) ' 

LEV - LEV+1 


C CHECK FOR ERROR. IF NOT ERROR. GO BACK AND DO THE NEXT LOOP 


I F ( ( V ( LEV ) - R ( LEV ) ) . G E . HMIN) GO TO 100 
225 IER - 129 

BB - SNGL(R(LEV) ) 

RETURN 

END 


SUBROUTINE CAUSS (M . C , RR) 

DIMENSION C(M.M) , RR(M) 

DD-0 . 0 
DO 23 L-l.M 
DO 21 K-l.M 

DD-DD+C ( L , K ) *C ( L , K ) 

CONTINUE 
: CONTINUE 

PD-SQRT(DD) 

DET-1.0 

MAXIMIZATION OF PIVOT ELEMENTS 

DO 101 L - 1 , M 
BIC-C(L.L) 

ABIG-ABS (C ( L , L) ) 

IBIG-L 

DO 15 I - L.M 
SIZE-C ( I , L) 

ASIZE-ABS ( C ( I . L) ; 

IF (ASIZE.LE.ABIG)GO TO 15 

BIG-SIZE 

IBIG-I 

15 CONTINUE 


DET-DET*BIG 

C SWAP ROUS FOR MAXIKUM OF C(L,L) 
IF (L.EQ.IBIG)GOTO 26 
DO 25 J-L.M 
CBIG-C ( IBIG , J ) 

C(IBIG, J)-C(L.J) 

C(L.J)-CBIC 

25 CONTINUE 
RBIG-RJU IBIG) 

RRdBIG)-RR(L) 

RR(L)-RBIG 

26 RR(L)-RR(L)/BIC 
DO 151 J-L.M 

151 C(L.J)-C(L.J)/BIG 
DO 201 K - l.M 
A-C(K.L) 

IF(K.EQ.L) GOTO 201 
DO 11 J-L.M 

11 C(K,J)-C(K,J)-C(L.J)*A 

RR(K)-RR(K) -RR(L)*A 

201 CONTINUE 

101 CONTINUE 

CTX-ABS (DET) /PD 

RETURN 

END 


APPENDIX C: FORTRAN PROGRAM FOR TVO DIMENSIONAL 
RADIATION/ CONDUCT I ON MODEL 

PARAMETER (ND - 100 . ND2 - 200 ) 

EXTERNAL AFUNC, THERM! . THERM 2 ADSIMP 
DIMENSION A ( ND . ND ) , B ( ND ) , T ( ND ) , IT ( 5*. 5 ) V(ND) 

S DI FF(ND) . AFT( 8 ) , ER( 2 ) , TC(ND) QT(ND) 

S WW f ND ) . EE ( ND ) . FT (f ND2 ) , D I S T ( ND ) 


SIGMA - 5 . 669E- 12 
TOL - L.OE-6 
IER - 0 


OPEN (6, FILE - 'CRUC1.DAT' 
INITIALIZE THE MATRICES 
DO 1 I - 1, nd 
BCI) - 0.0 


STATUS 


'NEW' ) 


TCI) - 0.0 
■-( I) - 0.0 
WW(I) - 0.0 
DO 1 J - 1 , ND 
A(I.J) - 0.0 
1 CONTINUE 

RANGE OF WAVE LENGTHS AND INITIAL GUESS TEMPERATURES 
WRITER.*) ’ENTER T(l) T(ND) ' 

READ(* , *) T( 1 ) , T(ND) 

WRITE(*,*) 'ENTER SL AND BL ' 

READ ( * , * ) SL, BL 
H - 1/ND 
DISK 1) - 0.0 
DO 6 I - 2, ND 

DIST(I) - DIST(I-l) + H 
CONTINUE 
COUNT - 0.0 

00 IF (COUNT. CT. 0.0) THEN 

IF ( FLX2 . GT . FLX1 ) THEN 
T(ND) - T(ND) + 3.0 
P - T(ND) 

ELSE 


T(ND) - T(ND) - 2.5 
P - T(ND) 

END IF 
ELSE 

WRITE ( 6 ,*) 'COUNTER - '.COUNT 
END IF 

DO 5 I - 1, 8 

AFT ( I ) - ADS IMP ( AFUNC , SL . BL . TOL . IER , I ) 

CONTINUE 

CALL RAD(ND ,T.B,H,SL,BL.ND, FLX1) 

WRITE ( * , *) 'FROM RADIATION FLUX AT JOINT - ' . FLXi 
DO 2 I - 1, ND 
FTC I) - T ( I ) 

CONTINUE 

CALL CONDUCT ( TT , T , H , HH , ND , FLXI , FLX2 ) 


102 


WRITE(* , *) 'FROM CONDUCTION FLUX AT JOINT - ' FLX2 
COUNT - COUNT + l 
IF(COUNT . GT . 100 . 0 ) GOTO 300 
C CONDITION FOR STOP 

I F( COUNT . GT . 0 . 0) THEN 
ER(2) - ER( 1) 

ER( 1 ) - ABS(FLX2 - FLX1 ) 

ERA - ABS ( ER( 1 ) • ER( 2 ) ) 

IF (ERA.GT. 1.0E-2) THEN 

WRITE(*,*) 'ERROR IN FLUXES '.ERA 
GOTO 200 
ELSE 

VRITE(*,*) * ERROR AND FLUX1 FINAL ' , ERA , FLX1 
END IF 
ELSE 
END IF 

300 WRITE(* , *) '# OF ITERATIONS - '.COUNT 

WRITE (*,*) 'ERROR IN FLUX AT JOINT - ' . ERA 
WRITE(*,*) ' FLUX AT RAD/COND JOINT - ' , FLX1 , FLX2 
DO 4 I - 1, ND2 

WRITE(*,*)I, ’ - ' . FT ( I ) 

4 CONTINUE 

DO 3 I - 1, 5 

DO 3 J - 1,5 

WRITE(* , *) I , ' , ' , J , ' - ' , TT( I , J ) 

3 CONTINUE 
STOP 
END 

C 

SUBROUTINE RAD ( ND , T , B , H , SL , BL , ND , FLX1 ) 

DIMENSION A(ND.ND) ,B(ND) ,T(ND) ,W(ND) .DIFF(ND) . AFT(8) 
$ TC(ND) ,QT(ND) ,UW(ND) , EE(ND) , FLXR(ND) , 

$ PT(ND) 

T(l) - Q 
T(ND) - P 
PT( I) - T(l) 

COUNT2 - 0.0 
SIGMA - 5 . 669E- 12 
TOL - 1.0E-6 
IER - 0 
WW(1) - 0.0 
WWW - 0.0 
DO 423 I - 2 , ND 
WW(I) - H + WWW 
WWW - WWW +H 
423 CONTINUE 

Z - (T(l) - T(ND) )/ND 
DO 4 I - 2, ND-1 
T( I ) - T(l) - 1*2 
PT( I ) - T(I) 

4 CONTINUE 

DO 41 I - 1 , ND 
QT(I) - T(I)/T(1) 

41 CONTINUE 



MM - 1 

C ADJUSTING GUESS VALUES 
IF(MM . GT . 0)GOTO 212 
37 DO 42 I - 2, ND* 1 

IF(PT(I) .GT.T(I) ) THEN 
PT ( I ) - T( I ) 

T(I) - T ( I ) + 0. 5*DIFF(I) 

QT( I ) - T(I)/T( 1) 

ELSEIF(PT(I) . LT.T(I) )THEN 
PT( I ) - T( I ) 

T(I) - T(I) - 0. 5*DIFF(I) 

QT(I) - T(I)/T(1) 

ELSE 
END IF 
42 CONTINUE 
212 Y-1.0 

DO 22 I - 1,ND 

TC(I) - 0 . 0 1*THERM1 ( T ( I ) ) 

22 CONTINUE 

C GENERATE MATRIX [A] and [B] 

A(l.l) - 1.0 

A(ND.ND-l) - TC(ND) + TC(ND-l) 

A(ND , ND) - -A(ND.ND-l) 

DO 2 I - 2, ND-1 

A(I.I-l) - (0.5*(TC(I-1) + TC(I) ) )/(H**2) 

A( I , I ) — (0.5*TC(I-1) + TC( I ) + 

$ 0.5*TC(I+1))/(H**2) 

A( I , 1+1 ) - (0 . 5*(TC ( I ) + TC(I+1)))/(H**2) 

2 CONTINUE 
C 

CALL FORCE (ND , QT , W , B , AFT , SL , BL , TOL , IER) 

B(l) - 1.0 

B(ND) - AFT ( 2 ) *SIGKA*(T(HD) *+4) 

CNTERM -(SIGMA*(T(1)**3))/AFT(3) 

DO 5 I - 2, ND-1 

B( I ) - - 1 . 0*(CNTERM*B( I ) ) 

5 CONTINUE 
C 

CALL GAUSS (ND , A, B) 

C 

DO 666 I - 1, ND 
QT(I) - B(I) 

T(I) - T(1)*B(I) 

666 CONTINUE 

DO 8 I - 1, ND 

W(I) - T(1)*B(I) 

WRITE(*,*)WW(I) , ' , ' ,W(I) 

8 CONTINUE 

FLXR(ND)- ( 3 . *T(ND) - 4.*T(ND-1) + T(ND-2) )/(2 .*H) 
FLX1 - 0 . 4*SIGMA*(T( 2 )**4 - T(l)**4) + FLXR(ND)* 

$ THERM1 (T(ND) ) 

DO 13 I- 2 , ND 

DIFF(I) - ABS (PT( I ) - T(I) ) 

13 CONTINUE 



CONDITION FOR STOP 

— WfVi ± ' 

COUNT 2 - COUNT 2 + L.O 
IF (COUNT2 . GT. 20 . 0) COTC ICC 
DO 14 I - 1. NO 

IF(DIFF(I) .GT.O. 5) GOTO 37 
CONTINUE 

i : : VRITEC*,*) ' total ITRS FOR RADIATION ’ . COUNT 2 

RETURN 

END 

SUBROUTINE CONDUCT (TT . T . H . HH , ND , FLX1 , FLX2 ) 
DIMENSION T(ND ) , U(ND) . DIFFC( 5 . 5) , BT( 25 ) 

$ BFLX( 2 ) ,TCC(5,5) , QT( 2 5 ) , UU( 25 ) , EE ( 2 5 ; . 

S FLXC( 25) . TT ( 5 , 5 ) . B B ( 2 5 ) , AA ( 2 5 . 2 5 ) 


GG - 1/5 
HH - 1/5 
DO 1 I - 1, 5 

TT(l.I) - T(ND) 

CONTINUE 
COUNT3 - 0.0 
SIGMA - 5.669E-12 
TOL - 1.0E-6 
IER - 0 
MM - 1 

C ADJUSTING GUESS VALUES 
I F ( MM . GT . 0 ) GOTO 212 

37 DO 42 N - 2, 5 

DO 42 M - 1. 5 

TT(N.M) - TT(N.M) • 0 . 1*DIFFC (N ,M)*T (1.1) 

ELSE 

TT(N,M) - TT ( N , M ) + 0 . 1*DIFFC( I )*TT( 1 . 1) 

END IF 

-2 CONTINUE 

212 Y-1.0 DO 22 N - 1,25 

DO 22 M - 1, 25 

TCC(N.M) - 0 . 01*THERM2 ( T ( N , M ) ) 

22 CONTINUE 

I GENERATE MATRIX [A] and IB] 

I R - GC**2/HH**2 

AB - ADSIMP ( AFUNC , SL , BL . TOL , IER , I ) 

DO 1 N - 2. 24 
DO 1 M - 2, 24 

A(N , M- 1) - ( 1/AB**3 ) *0 . 5* (TCC (N- 1 . M) +TCC ( N . M) ) 
A(N.M) - • ( ( AB**3 + HH*N ) * ( R/HH*N* ( AB** 3 ) ) ) * 

3 0. 5(TCC(N-1,M) - TCC(Nfl.M)) 

A ( N , M+l ) - (l/AB**3)*0.5*(TCC(Nfl.M)+TCC(N.M)) 
A(N-l.M) - (R/HH*N)*0. 5*(TCC(N-1.M)*TCC(N.M) ) 
A(N+1,M) - (R/HH*N)*0 . 5*(TCC(N+1 ,M)+TCC(N.M) ) 
CONTINUE 

C BOUNDARY CONDITIONS 

C 



DO 2 NN - 1, 5 

TT ( NN . 1 ) - T ( NN , 2 ) 

2 CONTINUE 

DO 3 NN - 1, 5 

TT ( NN , 5 ) - T(NN,4) 

2 CONTINUE 

DO 4 MM - I, 5 

TT(l.MM) - TTU.l) 

U CONTINUE 

DO 5 MM - 1, 5 

TT(5.MM) - TT(4,MM) 

5 CONTINUE 

r 

i - l 

DO 6 M - 1, 5 

DO 6 N - 1, 5 
BB( I ) - T(N , M) 

EE ( I ) - BB( I ) 

I - I + I 

6 CONTINUE 

C 

DO 55 I - I, 25 
BT( I ) - 0.0 
55 CONTINUE 

C 

DN - 25 

CALL GAUSS ( DN , A , BT) 

C 

C 

I - 1 

DO 66 M - 1, 5 

DO 66 N - I, 5 
T(N , M) - BT( I ) 

BB( I ) - BT( I ) 

I - I + 1 
66 CONTINUE 
C 

BFLX(2)-BFLX(1) 

FLXC(l)- ( - 3 . *BB( 1 ) + 4 . *BB ( 2 ) • BB ( 3 ) ) /( 2 . *HH) 
URITE(*,*) 'CONDUCTION flux at 1 - \FLXC(1) 
FLXC(25)- ( 3 . ★BB (25) - 4.*BB(24)+ 

$ BB(23))/(2.*HH) 

VRITE(*,*) 'CONDUCTION flux at 25 - \FLXC(25) 
FLX2 - FLXC(l) 

BFLX(l) - FLXC(25) 

FLX2 - THERM2 ( BB ( 1 ) ) *FLXC ( 1 ) 

DO 13 I- 2 , 25 

DIFFC(I) - ABS (BB ( I ) - EE(I)) 

13 CONTINUE 

C CONDITION FOR STOP 

WRITE(*,*) 'CONDUCTION Iterations - ' ,MM 
MM - MM+1 

COUNT 3 - COUNT 3 +1.0 


IF (COUNT 3 . GT . 20 . 0) GOTO 300 
DO 14 I - I, ND 

IF(DIFF(ND) .GT.1.0E-6) GOTO 37 
14 CONTINUE 

300 VRITE(* , *) ' TOTAL ITRS FOR CONDUCTION — ' COUNT 3 
RETURN 
END 

C 

SUBROUTINE FORCE (ND , ST , V , SS , AFT , SL , BL, TOL IER) 
DIMENSION ST(ND) , W(ND) , SS (ND) , AFT ( 8 ) 

REAL Y , XY , AFUNC 
C 

SIGMA - 5 . 669E- 12 
DO 1 I - 1, ND 
C 

SSE - 0.0 
C 

C GAUSSIAN QUADRATURE 
C 

DO 2 J - 1 , 2 
IF (J.EQ.l) THEN 

Y - -0.57735027 
LAMDA - SL 

ELSE 

Y - 0.57735027 
LAMDA - BL 

END IF 

C 

XY - 5 . 0E- 5*Y + 5.0E-4 
C 

A - AFT (1)*(ST(I) **4 ) 

B - ((1.0 - AFT(2) )*0 . 1281) / 

$ (2.0*0. 1281*(1.0 - AFT(2) ) 

$ + ( (AFT(l) )*SIGMA*(ST(I)**4) 

$ *( (ST(1)*1600)**4) ) ) 

C - AFT (3)*(ST(I) **4 ) 

D - ( AFUNC (1 .0,8) )*(ST(I)**4)* 

$ ( ( AFUNC (XY, 1) ) -AFT (4) ) 

F - 2 . 0*( 1 . 0 - (AFUNC (XY, 6))) 

G - ( AFUNC (XY ,6))*(ST(I) **4) 

H - 1.0/0.1281 
C 

A I LOG - 1.0 - XY 
SSE - SSE + 5.0E-5 * (4.0 * G 
$ - 2.0 * (AFUNC(XY, 2) ) * 

$ (G + 2 . 0*F*B* ( (-G - F*A) - H*(C-A) ) + A) 

$ + 4 . 0*B* ( - G - F*A - H*(C-A) ) 

$ * ( L.O + (0.5772 - 1.0 + alog(AlLOG) ) * 

$ (l.-XY) - ((1. -XY)**2)/2.0) - 2*D ) 

2 CONTINUE 

SS(I) - SSE 
1 CONTINUE 

RETURN 
END 


DO 666 I - 1, DN 
QT(I) - B(I) 

T ( I ) - B( I ) *T( 1) 

666 CONTINUE 

DO 8 I - 1, ND 

U(I) - T(1)*QT(I) 

WRITE(*,*)UU(I) , ' , \U(I) 

8 CONTINUE 
C 

REAL FUNCTION AFUNC (X, ID) 

REAL X 
INTEGER ID 

5 GOTO (10, 20, 30, <*0,50, 60, 70,80, 90) ID 

10 AFUNC - -0.5772 • alog(X) + X - (X**2)/4.0 
goto 90 

20 AFUNC “ 1.0 + (0.5772 • 1.0 + alog(X))*X - 
$ (X**2)/2 . 0 

goto 90 

30 AFUNC - 0.5 - X + 0 . 5* ( -0 . 5772 + 1.5 - alog(X) ) *X**2 
$ + (X**3)/6.0 

goto 90 

40 AFUNC - ((-0.0001 * (X**5 ) + 0 . 0031*(X**4) - 

$ 0 . 0523*(X**3) + 0.3834 * (X**2) - .9968 * X + 

$ .8177) * (1.003**4) * (-0.5772 - alog( 1 . 003E- 5 ) 

$ + 1.003E-5 - ( 1 . 003E- 5**2 ) /4 . 0) 

$ - ( -0 . 0001*(X**5) + 0.0031 * (X**4) - 0.0523 * 

$ (X**3 ) * 0.3834 * (X**2 ) - .9968 * x + .8177) * 

$ (1.003**4) *(-0.5772 - alog(X) + X - (X**2)/4.0)) 
goto 90 

50 AFUNC - (1.0 + (0.5772 - 1.0 + alog(X))*X - 

$ (X**2)/2.0) * ( -0 . 0001*(X**5) + 0.0031*(X**4) - 

$ 0 . 0523*(X**3)+ 0 . 3834*(X**2 ) - ,9968*X + .8177) 

goto 90 

60 AFUNC - -0.0001*(X**5) + 0 . 0031*(X**4) - 
$ 0 . 0523*(X**3) + 



$ 0 . 3834* (X**2 ) - 

. 9968*X + .8177 

70 

goto 90 

AFUNC - . 0004*X + .4432 


80 

goto 90 

AFUNC - 1. 0/(0. 5 - 1.0 + 

0 . 5*( -0 . 5772 + 1.5 


$ alog(l . 0) )*1 . 0**2 + (1, 

. 0**3)/6 .0) 

90 

RETURN 



END 

C 

C THERMAL CONDUCTIVITY OF ALUMINA 

C 

REAL FUNCTION THERMl(X) 

REAL X 

THERM1 - . 0019*X +6.0 

RETURN 

END 


C 

C 


o o o 


thermal conductivity of iron sample 


REAL FUNCTION THERM2 (X) 
REAL X 

THERM2 - . 0219*X +2.6 

RETURN 

END 


C- 


C 
C 
C 
C 
C 
C 
C 

C- - - 
C 


THIS FUNCTION CALCULATES THE INTEGRAL OF THE FUNCTION IN THE 
REAL FUNCTION 'AFUNC', USING ADAPTIVE QUADRATURE WITH SIMPSON'S 


A - LOWER BOUND OF THE INTEGRAL 
B - UPPPER BOUND OF THE INTEGRAL 
TOL - MAXIMUM ERROR OF THE QUADRATURE 

IER “ ER &OR. IER - 129 IF THE ALGORITHM COULD NOT CONVERGE. 


REAL FUNCTION ADS IMP (AFUNC , AA, BB , TOLL, IER.Q) 


REAL AA, BB, TOLL, AFUNC 
INTEGER Q 

DOUBLE PRECISION EST , R( 20) , V( 20) , FR( 20) , 

* FS(20) ,FT(20) ,FU(20) ,FV(20) ,ERR,S2F,SF, HMIN TOL 

LOGICAL FIN 

INTEGER LEV, IER 
C INITIALIZE 
IER - 0 
A - DBLE(AA) 

B - DBLE(BB) 

TOL - DBLE(TOLL) 

HMIN - (BB-AA)/2. DO** I 9 
EST - 0.D0 
R(l) - A 
V(l) - B 

IF (A.EQ.B) RETURN 
1 LEV - I 


FR(1) - DBLE (AFUNC ( SNCL( R(l) ) ,Q) ) 

FS(1) - DBLE(AFUNC ( SNGL( (V(l) *R(I) ) /4 . DO + R(I) 
FT(1) - DBLE(AFUNC ( SNGL( (V( I) -R(l) )/2 . DO + R(I) 
FU( 1 ) - DBLE(AFUNC( SNGL( (V(I) -R(l) )*. 75DO + R(l) 
FV(1) - DBLE (AFUNC (SNGL(V(1) ) ,Q) ) 

FIN - .FALSE. 


C 

C ENTER THE MAIN LOOP 
C 

100 
C 

C CALCULATE 
C 


IF (FIN) RETURN 


S2(F) AND SF( F) 


.Q) ) 
.Q> ) 
,Q) ) 


S2F - (V(LEV) -R(LEV) )/12 . D0*(FR(LEV)+FV(LEV)+2 . D0*FT(LEV) + 
4 . D0*( FS ( LEV) + FU ( LEV ) ) ) 

SF - (V(LEV) -R(LEV) )/6 . DO*(FR(LEV)+FV(LEV) + 4 . DO*FT(LEV) ) 



n o o 


* 4 ' J 


C GET THE TOLERANCE FOR THIS STAGE 

ERR - 2. DO*TOL* < V ( LEV 1 • R LEV \ ' / ( 3 . 4 

! !? THE ERROR IS TOO HIGH. CO TO THE NEXT SMALLER STACE 
C OTHERWISE. DO THIS STACE. 

IF(DABS(S2F-SF) . GT . ERR) GO TO 200 
EST - EST * S2F 
ADSIMP - SNCL(EST) 

C CHECK FOR ENDING CRITERIA. IF LEVEL - 1 AND V(LEV) - B. STOP 

r* 

IF ( ( LEV . GT . 1 ) . OR . ( V ( LEV ) . LT . B ) ) GO TO 150 
FIN - .TRUE. 

GO TO 100 
150 CONTINUE 

OTHERWISE, SET UP FOR THE NEXT STAGE 

R(LEV - 1) - V ( LEV ) 

FR(LEV • 1) - FV(LEV) 

FT (LEV - 1) - FU ( LEV - 1) 

FS ( LEV - 1 ) - DBLE(AFUNC(SNCL((V(LEV-1) -R(LEV-l) ) 

* /4 . DO+R ( LEV - 1 ) ) , Q ) ) 

FU(LEV-l) - DBLE ( AFUNC ( SNGL( ( V ( LEV - 1 ) -R( LEV - 1 ) )* . ? 5D0 - 

* R(LEV-l)) .Q) ) 
LEV - LEV - 1 
GO TO 100 

C IF IT DIDN’T WORK OUT, GO TO THE NEXT HIGHER LEVEL 
200 CONTINUE 


C SET UP THE BOUNDARIES AND FUNCTIONS FOR THE NEXT LEVEL CALCULATIONS 


★ 


★ 


IF (LEV. EQ. 20) GOTO 225 
R(LEV+1 ) - R(LEV) 

V(LEV+1) - (R(LEV)+V(LEV) )/2 . DO 
FR( LEV+1 ) - FR(LEV) 

FT ( LEV+1 ) - FS (LEV) 

FV( LEV+1) - FT (LEV) 

FS ( LEV+1 ) - DBLE (AFUNC ( SNGL( (V( LEV+1) -R( LEV+1) ) /4 . DO 

R( LEV+1) 

FU ( LEV+1 ) - D BLE( AFUNC ( SNGL( (V( LEV+1) *R( LEV+1) )* . 75DO 


LEV - LEV+1 


R( LEV+1; 


Q) .> 
5 ) ) 


C CHECK FOR ERROR. IF NOT ERROR, GO BACK AND DO THE NEXT LOOP 


IF( (V( LEV) -R( LEV) ) . GE . HMIN) GO TO 100 
225 IER - 129 

BB - SNGL(R(LEV) ) 

RETURN 

END 


SUBROUTINE GAUSS (M , C , RR) 

DIMENSION C(M,M) ,RR(M) 

DD-0 . 0 
DO 23 L-l.M 
DO 21 K-l.M 

DD-DD+C (L , K) *C (L , K) 

21 CONTINUE 
23 CONTINUE 

PD-SQRT(DD) 

DET-1.0 

C MAXIMIZATION OF PIVOT ELEMENTS 

DO 101 L - 1 , M 
BIG-C ( L , L) 

ABIG-ABS(C(L, L) ) 

IBIG-L 

DO 15 I - L.M 
SIZE-C ( I , L) 

ASIZE-ABS (C(I ,L) ) 

IF (ASIZE.LE.ABIG)GO TO 15 

BIG-SIZE 

IBIG-I 

15 CONTINUE 
DET-DET*BIG 

C SWAP ROWS FOR MAXIMUM OF C(L,L) 

IF (L.EQ.IBIG)GOTO 26 
DO 25 J-L.M 
CBIG-C( IBIG , J ) 

C(IBIG, J)-C(L, J) 

C(L,J)-CBIG 

25 CONTINUE 
RBIG-RR(IBIG) 

RR(IBIG)-RR(L) 

RR(L)-RBIG 

26 RR(L)-RR(L)/BIG 
DO 151 J-L.M 

151 C(L,J)-C(L,J)/BIG 
DO 201 K - 1 ,M 
A-C (K, L) 

IF(K.EQ.L) GOTO 201 
DO 11 J-L.M 

11 C(K,J)-C(K,J)-C(L,J)*A 

RR(K)-RR(K)-RR(L)*A 
201 CONTINUE 
101 CONTINUE 

CTX-ABS (DET)/PD 

RETURN 

END 
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